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1.  SUMMARY 


This  report  describes  work  carried  out  under  the  Air  Force  Research  Laboratory's 
basic  research  task  in  optical  remote-sensing  signatures,  entitled  Optical  /  Infrared 
Signatures  for  Space-Based  Remote  Sensing.  This  work  was  carried  out  during  the 
period  20  Mar  2000  -  30  Sep  2007.  It  was  supported  by  funding  from  AFOSR's  Physics 
and  Electronics  Directorate  (AFOSR/NE),  in  the  portfolio  of  program  manager  Kent 
Miller. 

The  objective  of  the  program  was  to  increase  understanding  of  remotely  sensed 
signatures,  particularly  as  seen  by  spacebome  optical/infrared  (IR)  sensors.  Initially,  the 
emphasis  was  on  signatures  of  importance  to  missile  defense,  while  later  studies  focused 
on  signatures  relevant  for  sensing  resident  space  objects  (RSOs)  and  related  to  space- 
weather  applications.  Much  of  the  work  focused  on  the  environment  of  space  objects  and 
missiles,  especially  that  constituting  the  mesosphere  and  lower  thermosphere  (MLT). 

Much  of  the  work  carried  out  was  based  on  two  noteworthy  satellite  sensor  databases. 
These  are  the  unique  but  episodic  data  from  the  DoD  MSX  (Midcourse  Space 
Experiment)  satellite  and  the  extensive  global  data  from  the  NASA  TIMED  satellite. 

Both  satellites  carried  advanced  IR  sensors  in  addition  to  other  ultraviolet  (UV)  and  other 
sensors.  In  the  case  of  MSX,  this  was  the  SPIRIT  III  spectrometer  and  multichannel 
radiometer  array,  while  TIMED  bore  the  SABER  multichannel  radiometer  into  orbit. 
MSX,  launched  in  1996,  has  proven  to  be  a  rich  source  of  information  on  a  number  of 
phenomena  of  importance,  with  the  emphasis  on  atmospheric  gravity  waves  (GWs)  and 
on  auroras  during  this  reporting  period.  TIMED  was  launched  in  the  second  year  of  this 
effort,  so  the  period  of  performance  coincides  quite  nicely  with  the  first  six  years  of 
SABER  data.  TIMED/SABER  research  focused  on  (1)  developing  and  validating 
retrieval  algorithms,  especially  algorithms  for  retrieving  temperature,  pressure,  and  CO2 
volume  mixing  ratio  (vmr),  and  (2)  conducting  science  investigations  based  on  the 
SABER  data.  The  science  investigations  have  included  studies  of  temperature-inversion 
(TI)  layers,  tides,  radiance  layers,  GWs,  solar  storms  and  auroras,  and  OH  interannual 
variability. 

The  structure  of  the  radiance  environment  is  particularly  important  when  one  is 
seeking  a  space  object,  since  it  can  reduce  the  detection  probability  and/or  increase  the 
probability  of  false  alarm.  Radiance  structures  can  often  be  explained  by  understanding 
the  structure  of  the  underlying  medium.  A  special  emphasis  was  placed  in  this  period  on 
understanding  localized  structures  of  the  MLT.  These  include  vertically  localized 
structures  such  as  TI  layers  and  horizontally  localized  structures  such  as  fronts.  One 
particularly  important  type  of  MLT  front  is  the  mesospheric  bore,  identified  earlier  by 
AFRL  team  members  and  studied  further  during  the  present  period. 

Progress  was  also  made  in  the  current  period  on  the  better  studied  forms  of  MLT 
structure,  including  atmospheric  GWs,  tides,  and  turbulence.  One  investigation  observed 
and  characterized  concentric  GWs  and  upward  lightning  from  cloud  tops  known  as 
sprites  generated  simultaneously  by  the  same  thunderstorm. 

The  ability  to  predict  the  radiative  environment  depends  on  a  deep  knowledge  of 
radiative  transfer.  In  the  MLT  region  the  atmosphere  is  never  in  local  thermodynamic 
equilibrium  (LTE).  So  radiative  transfer  becomes  much  more  complicated,  and  non-LTE 
methods  are  required.  Progress  on  non-LTE  radiative-transfer  methods  occurred  during 
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this  period,  especially  methods  for  radiative  transfer  in  structured  atmospheres.  This 
knowledge  was  applied  to  improving  remote-sensing  retrievals.  Team  members  also 
collaborated  as  critics  and  reviewers  with  the  authors  of  the  first  major  book  to  be 
published  on  non-LTE  atmospheric  radiative  transfer. 

During  this  reporting  period  team  members  have  participated  in  organizing  an 
ongoing  remote-sensing  meeting  and  giving  it  a  new  focus  on  remote  sensing  of  the 
middle  and  upper  atmosphere.  They  also  applied  their  knowledge  of  radiative  transfer  to 
coupling  atmospheric  and  foliage-canopy  radiative-transfer  models  in  order  to  investigate 
possibilities  for  passive  foliage  penetration  from  a  satellite  sensor. 

Toward  the  end  of  the  period,  studies  began  on  exploiting  ambient  earth-atmosphere 
radiation  to  illuminate  RSOs  and  on  various  methods  to  enhance  detection  and 
characterization  of  space  objects. 

The  team  has  collaborated  with  a  number  of  programs  including  the  Columbia 
Space  Shuttle  Accident  Investigation,  the  HAA  program  and  the  Near-Space  Study, 
C/NOFS,  and  Maui  MALT. 

The  work  of  the  team  was  recognized  by  the  granting  of  the  Harold  Brown  Award  to  a 
team  member  by  the  Secretary  of  the  Air  Force  and  by  the  conferring  of  the  Guenther 
Loeser  Award  and  the  Marcus  O'Day  Award  to  team  members  by  the  Battlespace 
Environment  Division  (AFRL/RVB). 
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2.  INTRODUCTION 


Much  of  the  work  in  this  task  was  focused  on  the  environment  of  space  objects  and 
missiles,  especially  the  MLT  environment,  with  a  special  interest  in  a  highly  structured 
environment.  The  reason  for  this  focus  was  that  a  structured  environment  can  mimic  a 
target  of  interest  and  be  confused  with  it,  resulting  in  a  false  alarm,  or  can  act  as 
unintentional  camouflage  enabling  a  target  to  hide,  resulting  in  a  missed  detection. 

The  MLT  region  is  at  once  the  upper  limit  of  the  atmosphere  and  the  lowermost 
region  of  space.  As  such,  it  is  a  transition  region  where  a  wholesale  change  of  regimes 
takes  place,  including  from  turbulent  to  free  molecular  flow  and  from  LTE  to  non-LTE 
radiative  transfer.  Energy  is  deposited  by  energetic  solar  radiation  and  by  precipitating 
electrons  and  ions.  Photochemical  reactions  and  photochemical  heating  start  to  become 
very  important,  and  the  ionosphere  starts  to  develop.  As  a  consequence,  free  radicals  and 
ions/electrons  become  significant  species. 

At  the  same  time  the  region  is  notoriously  difficult  to  study,  both  experimentally  and 
theoretically.  Since  it  is  above  the  range  of  high-altitude  balloons  ands  aircraft  and  below 
the  altitude  where  stable  satellite  orbits  are  possible,  in-situ  measurements  can  rarely  be 
made  there,  and  most  of  our  knowledge  of  the  region  must  come  from  remote  sensing. 
Moreover,  models  are  difficult  because,  as  discussed  above,  virtually  everything  is 
important  there  and  no  effects  can  be  neglected. 

In  order  to  address  these  issues  we  have  developed  an  approach  emphasizing 
collaboration  with  field  and  space  experiments,  data  analysis,  modeling,  and  theory.  The 
approach  was  multi-pronged  and  involved  a  number  of  generic  tasks,  including  the 
following: 

•  Investigate,  specify,  and  exploit  optical  /  1R  spatiotemporal  and  spectral 
signatures  for  space-based  remote  sensing  of  space  objects  and  their  environment 

•  Construct  first-principles  models  of  phenomena  contributing  to  important  optical  / 
1R  signatures  measured  by  existing  instruments 

•  Extrapolate  measurements  using  models  from  the  measured  parameter  subspace 
to  the  total  range  of  parameters  of  AF  interest 

•  Validate  models  of  optical  /  1R  signatures  by  participation  in  satellite  and  ground- 
based  field  measurements 

•  Search  for  new  phenomena  to  exploit  and/or  mitigate  and  new  techniques 
applicable  to  increasing  awareness  of  space  objects  and  their  environment. 

We  have  expended  much  effort  to  analyze,  understand,  and  model  data  from  the 
pioneering  MSX  and  TIMED  space  experiments,  and  have  collaborated  as  well  with 
ground-based  observers  and  other  analysts  and  modelers. 

The  results  have  been  notable  advances  in  knowledge  of  a  number  of  aspects  of 
remote-sensing  signatures.  We  review  these  advances  in  the  following  sections.  Sections 
3  and  4  deal  with  the  MSX  and  TIMED/SABER  results,  respectively.  In  Section  5  we 
discuss  research  in  radiance  structure  and  how  it  is  related  to  the  structure  of  the 
underlying  medium,  concentrating  especially  on  localized  structure.  Research  in 
radiative  transfer  and  applications  to  remote  sensing  are  discussed  in  Section  6.  Section  7 
treats  preliminary  research  in  surveillance  of  space  objects.  In  Sections  8  and  9  we 
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discuss  collaborations  and  conclusions.  The  report  ends  with  references,  a  list  of  refereed 
and  unrefereed  publications  under  the  task,  and  a  list  of  acronyms  and  symbols.  Further 
details  of  closely  related  work  were  described  in  *Donatelli  et  al.  [2005], 


3.  MIDCOURSE  SPACE  EXPERIMENT  (MSX) 

3.1.  MSX  Mission  and  Data 

The  Midcourse  Space  Experiment  (MSX)  program  was  planned  and  implemented  to 
support  the  design  requirements  for  the  next  generation  of  space-based  IR  surveillance 
systems  and  to  provide  functional  demonstrations  of  target  detection  against  terrestrial, 
atmospheric,  and  celestial  backgrounds.  The  satellite  was  launched  on  24  April  1996 
from  Vandenberg  AFB  on  a  Delta  II  rocket  into  a  nearly  sun-synchronous  orbit  at  an 
altitude  of  900  km.  As  part  of  an  extensive  series  of  measurements,  MSX  successfully 
completed  150  episodic  measurements,  each  nominally  20  to  30  minutes  in  duration,  of 
IR  atmospheric  radiance  and  spatial  structure  in  earthlimb  scenes  at  tangent  heights 
ranging  from  30  to  150  km.  The  principal  sensor  onboard  was  the  SPIRIT-III  (Spatial 
Infrared  Imaging  Telescope  III)  radiometer  and  Michelson  interferometer/spectrometer, 
developed  by  Utah  State  University's  Space  Dynamics  Laboratory,  which  made 
measurements  in  a  number  of  mid-wave  IR  (MWIR)  and  long-wave  IR  (LWIR)  bands. 
The  sensor  suite  also  included  the  UVISI  SPIMs  (spectral  imagers)  and  narrow- field  and 
wide-field  imagers,  developed  by  Johns  Hopkins  University  Applied  Physics  Laboratory 
and  measuring  from  the  visible  into  the  ultraviolet  (UV).  IR  measurements  continued 
until  the  cryogenic  supply  of  solid  hydrogen  was  depleted  on  25  February  1997.  In  one 
very  productive  type  of  experiment,  high-resolution  infrared  images  were  recorded  in  a 
pushbroom  scan  mode  by  maintaining  the  SPIRIT  III  line-of-sight  (LOS)  at  90°  to  the 
LOS  in  the  plane  normal  to  the  orbital  plane.  The  MSX  program  objectives,  sensors,  and 
experiments  are  described  further  by  Mill  et  al.  [1994]  and  O  'Neil  et  al.  [1994], 


3.2.  Atmospheric  Gravity  Waves 

Efforts  under  this  work  unit  dealing  with  GWs  build  on  earlier  pioneering  papers  by 
team  members  on  observing  GWs  seen  in  the  IR  from  earth  orbit  and  modeling  the 
process  [Dewan  et  al.,  1998;  Picard  et  al.,  1998],  The  first  paper  analyzes  two  examples 
of  circular  concentric  GWs  seen  in  MSX  Bl-band  (4.22-4.36  pm)  below-the -horizon 
(BTH)  near-limb  scenes  and  shows  that  the  waves  are  radiated  from  point  sources 
associated  with  compact  convective  storms.  The  second  paper  calculates  the  spectrally 
resolved  contribution  functions  (CFs)  for  C02(v3)  emission  at  4.3  pm  for  various  LOS, 
including  ATH  (above-the-horizon)  and  various  BTH  cases,  both  nadir  and  near-limb 
(sub-limb).  By  integrating  over  the  B1  filter  function,  the  paper  shows  that  the  MSX  Bl- 
band  CF  peaks  near  45  km  altitude  and  has  no  contribution  from  below  the  upper 
stratosphere.  These  and  other  MSX  results  on  GWs  were  reviewed  in  O'Neil  et  al. 

[1998], 
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Large  spatial  fluctuations  and  structure  have  been  found  in  MSX  CO2  1 5 -pm  earthlimb 
data  at  tangent  heights  of  100-115  km  (Figure  1).  These  fluctuations  are  not  present  in 
other  bands  such  as  CO2  4. 3 -pm  and  have  important  implications  for  LWIR  backgrounds 
as  well  as  for  the  lower-thermospheric  energy  budget.  The  fluctuations  are  most  likely 
due  mostly  to  natural  variations  in  the  radiance  layers,  or  "knees",  occurring  near  110-115 
km  (see  Section  4.5),  coupled  with  the  passage  of  GWs.  The  implications  for 
atmospheric  cooling  rates  are  especially  significant  (Figure  2),  with  substantial 
differences  of  up  to  50K/d  in  the  cooling  rates  associated  with  the  three  colored  profiles 
of  Figure  1  in  the  100-km  altitude  region  [*  Picard  et  al.,  2002] 1 .  Radiance  knees  are 
discussed  further  in  Sections  4.5  and  5.1.2. 


EL0201000134  Band  D  Data 
30-sec  Intervals,  UT  14:04-14:28 


Figure  1.  MSX  limb-radiance  profiles  taken  30  s  apart  in  a  view  direction  orthogonal  to  the 
velocity  vector,  showing  considerable  structure  due  to  radiance  layers  in  the  90-110  km  region. 

The  observed  radiance  from  the  MSX  SPIRIT-3  radiometer  contains  spatial  structure 
down  to  the  scale  of  hundreds  of  meters.  This  spatial  structure  results  from  local 
fluctuations  in  the  temperature  and  densities  of  the  radiating  states  of  the  emitting 
molecular  species,  as  well  as  fluctuations  in  radiation  transport  from  the  emitting  regions 
to  the  observer.  A  portion  of  this  database  has  been  analyzed  to  obtain  statistical 
parameters  characterizing  stochastic  spatial  structure  in  the  observed  radiance  [*. Brown  et 
al,  2002],  The  analysis  focused  on  five  scenes  in  the  CO2  4.3-pm  radiating  region  (4.22- 
4.36  pm,  MSX  B1  band)  and  one  scene  in  the  atmospheric  window  region  (6.8-10.8  pm, 
MSX  A  band).  The  scenes  included  both  ATH  (limb)  scenes  and  various  BTH  scenes, 


1  An  asterisk  (*)  indicates  an  unrefereed  publication  by  team  members  in  the  period  of  performance  of 
this  effort  (2000-2007),  while  a  double  asterisk  (**)  indicates  a  refereed  publication.  The  asterisks  also 
label  these  publications  in  the  list  of  References.  Team-member  publications  are  summarized  in  the 
Publications  section,  with  subsections  for  **Refereed  and  *Unrefereed  publications. 
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both  nadir  and  near-limb  (sub-limb).  Using  simple  models,  the  observed  statistics  were 
shown  to  agree  with  prior  observations  and  theoretical  models  of  stochastic  spatial 
structure  generated  by  superposed  and  interacting  GWs  for  special  viewing  geometries. 
The  SHARC  (Strategic  High  Altitude  Radiance  Code)  model  [ Sundberg  et  al,  1995]  had 
been  extended  (SHARC  4)  [Gruninger  et  al.,  1995,  1998,  2004]  to  predict  the  statistics  of 
stochastic  fluctuations  in  infrared  radiance  from  the  statistics  characterizing  temperature 
fluctuations  in  the  middle  and  upper  atmosphere  for  arbitrary  viewing  geometries. 

SHARC  model  predictions  were  compared  with  MSX  data  and  shown  to  be  in  generally 
good  agreement.  Measured  radiance  power  spectra  were  characterized  by  power  laws 
over  certain  ranges  of  wavenumber  which  varied  in  a  systematic  way  with  nadir  angle. 


C02  15  Micron  Cooling  Rate 


Figure  2.  C02 15-pm  cooling  rates  calculated  from  the  three  correspondingly  colored  scans  of 
Figure  1  with  the  help  of  runs  of  the  TIME-GCM  model  [Roble  and  Ridley,  1994],  showing  the  large 
variation  in  cooling  rate  due  to  the  presence  of  radiance  layers. 

Some  of  the  IR  radiance  power  spectra  showed  breakpoints  where  changes  of  power-law 
slope  occurred.  Both  geophysical  and  observational  explanations  for  breakpoints  in  the 
spectra  originating  from  an  atmosphere  perturbed  by  superposed  GWs  were  proposed 
[ *Picard  et  al.,  2003a].  Previous  works  on  the  spectra  of  known  sources  of  GWs  and 
wave-filtering  mechanisms  by  M.  J.  Alexander  and  others  [Alexander,  1998;  Pandya  and 
Alexander,  1999;  Piani  et  al.,  2000]  and  on  the  general  expected  characteristics  of  GW 
power  spectra  [Dewan,  1979,  1997;  Cho  and  Lindborg,  2000]  were  reviewed.  Then  the 
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power  spectra  of  LOS  radiance  fluctuations  emitted  and  self-absorbed  by  an  atmosphere 
perturbed  by  GWs  were  considered,  including  the  shape  of  the  spectrum  and  the  spectral 
slopes.  Examples  of  radiance  spectra  from  gravity-wave-perturbed  atmospheres  that 
have  two  different  slopes,  with  a  steeper  slope  at  large  wavenumber,  were  analyzed,  and 
possible  mechanisms  were  discussed.  The  effect  of  latitude  and  season  on  the  4. 3 -pm 
fluctuations  was  also  considered. 

Although  it  is  clear  that  explanation  of  radiance  power-spectral  breakpoints  requires 
both  geophysical  and  observational  input,  future  research  will  be  needed  to  complete 
understanding  of  their  nature  and  properties.  More  emphasis  will  need  to  be  placed  on 
understanding  spectra  of  orographically  generated  waves  (mountain  waves)  and,  in 
addition,  the  investigation  will  need  to  understand  the  behavior  of  higher  moments, 
probability  densities  and  exceedances. 


3.3.  Auroras 

The  MSX  SPIRIT-3  radiometer  window  band,  band  A  (6.8-10.8  pm)  includes 
emission  from  ro vibrational  bands  of  CO2,  H2O,  and  O3.  Under  quiescent  conditions,  it 
does  not  include  any  emission  from  NO.  The  lower-thermospheric  NO  quiescent 
rovibrational  airglow  emission,  largely  excited  by  O  collisions,  occurs  almost  entirely  in 
the  first  fundamental  band  (transition  v-1—  >v"=0).  However,  auroral  electrons  induce 
chemiluminescent  reactions  between  N(2D)  or  N(4S)  and  O2,  leading  to  production  of 
higher  states,  v'>l,  in  the  fundamental  band  whose  emission  is  transmitted  in  part  by  the 
band-A  filter.  On  10  Nov  1996,  MSX  saw  enhanced  emission  in  the  A  band  at  night 
between  110  and  130  km  tangent  height,  and  it  was  shown  that  the  enhanced  emission 
resulted  from  NO(v)  produced  by  strong  auroral  dosing  [**Sharma  et  al.,  2001],  Using 
the  auroral  electron  energy  flux,  the  location  of  the  dosed  region,  and  the  MSX  LOS, 
along  with  the  N2+  First-Negative  band  emission  at  391 .4  nm  measured  onboard  by  the 
UVISI  SPIMs,  the  3-D  energy  deposition  rate  was  derived.  This  rate  was  input  to  the 
SHARC  auroral  model  [Sundberg  et  al.,  1995],  which  output  radiance  results  in  good 
agreement  with  MSX  measurements. 

A  unique  detailed  study  [** O’Neil  et  al.,  2007]  resulted  from  an  auroral  event 
observed  on  9  Feb  1997  by  MSX.  The  SPIRIT-3  radiometer  on  MSX  observed  enhanced 
4.3-pm  emission  from  a  very  well  characterized  aurora  over  the  Barents  Sea  near 
19:30UT,  in  conjunction  with  observations  by  the  Polar  and  FAST  satellites,  the 
Loparskaya  ground  site,  and  the  UVISI  ultraviolet  and  visible  spectrometers  aboard 
MSX.  The  MSX  track  and  the  locations  of  the  coordinated  measurements  are  shown  in 
Figure  3.  Measurements  of  the  auroral  location,  form,  spatial  extent,  and  dosing 
conditions  were  used  to  specify  the  component  of  auroral  4.3-pm  radiance  due  to  the 
time-delayed  and  optically  thick  C02  v3  (00°l-00°0)  radiation.  An  analysis  based  on  the 
AARC  (Auroral  Atmospheric  Radiance  Code)  IR  auroral  radiance  model  [Winick  et  al., 
1987]  indicates:  (1)  the  emission  originates  from  near  to  and  beyond  the  tangent  point; 

(2)  the  optically  thick  C02  v3  radiation  is  largely  self-absorbed  by  the  intervening 

atmosphere;  and  (3)  the  auroral  enhancement  is  predominantly  due  to  NO+  Av  =  1 
rovibrational  transitions.  In  addition,  the  analysis  indicates  that  the  previously  reported 
laboratory  result  for  the  NO+  v  >  1  vibrational  yield  from  the  reaction,  N+  +  02  (reaction 
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1)  is  insufficient  to  account  for  the  observed  4.3-pm  emission.  In  order  to  explain  the 
current  results,  we  proposed  that  there  is  additional  production  from  the  reaction,  N2+  + 

O  (reaction  2),  forming  NO+(v)  in  vibrational  levels  v  =  0,  1,  and  2  with  relative 
populations  of  approximately  0.25,  0.5  and  0.25,  respectively.  The  AARC  simulations  for 
CO2  alone  and  in  combination  with  an  NO+  contribution  are  shown  in  Figure  4.  The 
combined  production  processes  yield  an  energetic  electron- induced  radiant  efficiency  of 
0.56  ±0.18  photons  per  auroral  ion-pair  for  NO+  Av  =  1  emission  at  altitudes  equal  to  or 
greater  than  112  km.  This  work  enables  the  reliability  of  AF  predictive  models  of  auroral 
infrared  radiance  to  be  improved.  In  addition,  it  improves  our  knowledge  of  the  auroral 
energy  budget,  where  4.3-pm  emission  plays  a  major  role. 


Figure  3.  The  northern  polar  region  and  locations  of  the  MSX  satellite  and  coordinating  sensors  on  9  February  1997 
near  19:30  UT.  The  day-night  terminator  at  19:27  UT  is  indicated  for  altitudes  of  0,  95  and  900  km  and  the  orbital 
tracks  of  the  MSX  (red  trace)  and  FAST  (white)  satellites  are  shown  at  selected  times  in  the  period  from  19:25  to 
19:48  UT.  The  end  point  of  the  LOS  (yellow)  of  the  MSX  sensors  is  the  location  of  the  bore-sight  tangent  point  at  an 
altitude  of  95  km.  The  nadir  point  of  the  Polar  spacecraft  and  the  75°  half-angle  field-of-view  (FOV)  of  the  O  ( 1 1)) 
all-sky  camera  at  Loparskaya,  projected  to  250  km,  are  also  shown.  Beside  the  map  are  two  images  of  the  aurora 
from  the  MSX  UVISI  ultraviolet  sensors  at  19:29:14UT.  The  top  image  is  from  the  narrow-FOV  imager,  while  the 
bottom  one  is  from  the  wide-FOV  imager. 
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Rodtanca  Qmtta  on4  r1) 

Figure  4.  This  panel  shows  the  total  AARC  model  radiance  profile  (red  curve)  for  the  5-minute 
dosing  case  and  its  NO+  (blue  curve)  and  CO2  (green  curve)  components.  The  CO2  radiance 

profiles  for  dosing  times  of  1  and  3  minutes  are  also  shown.  The  NO+  emission  is  prompt  and, 
hence,  independent  of  the  dosing  history.  The  black  curve  is  the  MSX  B2-band  radiance 
profile. 


3.4.  Atmospheric  Cooling  Rates 

Atmospheric  radiative  cooling  rates,  important  for  the  atmospheric  energy  budget  and 
dynamic  models,  have  been  inferred  from  MSX  data  [** Picard  et  al.,  2002],  CO2  15-pm 
emission  is  primarily  responsible  for  cooling  the  atmosphere  from  the  tropopause  up  to 
about  120  km.  The  cooling  rates  can  become  very  large  above  90  km,  as  shown  in 
Figure  2.  In  addition,  they  are  quite  variable,  depending  a  great  deal  on  the  structure, 
including  the  fine  structure  of  the  atmospheric  temperature  and  species  profiles. 


4.  TIMED/SABER 

4.1.  TIMED/SABER  Mission  and  AF  Contribution 

The  TIMED  satellite  [Yee  et  al.,  2003]  is  the  first  remote-sensing  satellite  ever 
designed  to  explore  systematically  the  “ignorosphere”,  the  region  between  60  and  1 80 
km,  above  the  range  of  in  situ  measurements  by  high-altitude  balloons  and  below  the 
range  of  in  situ  measurements  by  satellites  in  high-eccentricity  orbits.  SABER  (Figure  5) 
[Russell  et  al.,  1999]  is  the  main  instrument  (one  of  four)  on  TIMED  and  consists  of  a 
100%-duty-cycle,  ten-channel  IR  limb-scanning  radiometer,  measuring  radiance  bands 
from  1.27  pm  [(^('Ag,  0-0)]  to  15  pm  [CC>2(v2)]  as  a  function  of  tangent  height  (Level  1 
products).  These  radiance  channels  are  inverted  to  yield  retrieved  atmospheric 
temperature,  pressure,  minor-species  profiles,  and  radiative  and  chemical  heating/cooling 
rate  profiles  (Level  2  products)  as  a  function  of  altitude,  as  well  as  further  more  highly 
processed  products.  Typically,  the  altitude  range  is  from  20  to  180  km  or  so,  where  the 
signal  level  reaches  the  noise.  In  some  cases,  such  as  NO  under  auroral  conditions,  the 
radiance  data  extends  to  considerably  higher  altitude  (up  to  300  km). 
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Figure  5.  TIMED  satellite  showing  position  of  four  instruments  including  SABER. 


SABER  logo  indicates  coonerating  entities  including  AFRL. 


TIMED  was  launched  on  7  Dec  2001  from  Vandenberg  AFB  into  a  near-circular  600- 
km  orbit  after  many  years  of  planning  and  of  launch  delays.  It  has  been  taking  data  since 
Jan  2002  and  is  now  completing  its  sixth  year  of  successful  operation.  The  TIMED 
mission  is  managed  by  NASA  Goddard  Spaceflight  Center  and  Johns  Hopkins  Applied 
Physics  Laboratory,  while  NASA  Langley  Research  Center  and  Hampton  University 
manage  SABER.  The  SABER  instrument  has  operated  with  a  100%  duty  cycle  since  the 
inception  of  the  mission  and  is  still  performing  near- flawlessly.  It  was  designed  and 
fabricated  by  the  Space  Dynamics  Laboratory  of  Utah  State  University. 

Every  one  of  the  SABER  channels  is  in  the  non-LTE  regime  over  at  least  part  of  its 
altitude  range.  This  is  true  for  both  the  target  species  and  some  of  the  interfering  species. 
Consequently,  the  retrieval  process  was  a  daunting  task  that  had  never  before  been 
attempted  on  this  scale.  An  essential  component  of  the  retrieval  codes  is  the  non-LTE 
forward  model  that  predicts  the  radiance  from  a  specified  atmosphere.  In  1995  AFRL 
was  asked  to  join  the  SABER  Science  Team  because  of  its  world-class  expertise  in  non- 
LTE  radiative  transfer.  AFRL  personnel  are  regular  participants  in  meetings  of  the 
SABER  Science  Team  and  also  of  the  parent  TIMED  Science  Working  Group.  See,  for 
example,  *  Picard  et  al.  [2004], 

Of  particular  relevance  was  AFRL's  longtime  experience  with  CO2  non-LTE 
emissions.  These  include  two  SABER  target  bands— the  V2  band  at  15  pm  and  the  V3 
band  at  4.3  pm — as  well  as  the  CO2  laser  bands  at  9.4  pm,  which  are  an  important 
interfering  emission  in  the  SABER  O3  channel  at  9.6  pm.  AFRL  team  members  have 
helped  to  design,  test,  and  validate  the  innovative  SABER  retrieval  algorithms  (for 
example,  **Mertens  et  al.  [2001]).  SABER  data  is  validated  by  comparison  with  models 
and  with  ground-based,  airborne,  and  other  space-based  data.  The  latter  include  data 
from  the  old  reliable  NASA  UARS/HALOE  probe,  from  MSX,  and  from  the  new 
ENVISAT-1  /  MIPAS  probe  launched  by  ESA  in  2003.  This  AFRL  involvement  ensures 
the  quality  of  the  unique  SABER  data  set  and  its  availability  and  utility  for  AF  purposes. 
It  benefits  especially  AF  models  and  codes  of  upper-atmospheric  and  space  radiative 
structure.  In  addition,  AFRL  team  members  analyze  data  for  significant  effects, 
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dependencies,  and  trends.  The  collects  to  date  have  yielded  remarkable  data,  increasing 
AF  knowledge  of  the  high-altitude  radiative  environment  and  of  the  underlying 
background  atmosphere  and  its  structure  manyfold.  The  database  has  gone  through 
several  versions,  the  current  version  being  vl.06  with  vl.07  soon  to  be  released.  The 
availability  of  global  SABER  retrievals  of  kinetic  temperature  (Tk),  pressure,  and  mixing  ratios 
of  C02  as  well  as  03  and  H20  on  a  24/7  basis  should  improve  model  atmospheres  such  as  SAG 
(SHARC  Atmosphere  Generator)  [ Adler-Golden ,  1993]  for  DoD  radiance  codes.  For  the  first 
time,  SABER  data  organized  by  latitude,  month,  and  time  of  day  will  provide  a  capability  to  the 
codes  to  characterize  intrinsic  variability.  Specific  examples  of  investigations  involving  AFRL 
personnel  are  in  the  following  paragraphs. 


4.2.  Temperature  and  CO2  Mixing  Ratio  Retrieval 

The  SABER  temperature  retrieval  is  critical,  since  most  of  the  other  retrievals  depend 
on  knowing  the  kinetic  temperature,  Tk,  well.  Also,  the  retrieval  over  the  target  altitude 
range  (20-130  km)  is  a  challenging  retrieval  that  is  in  the  non-LTE  regime  in  the 
mesosphere  and  above  and  had  never  been  attempted  before.  The  retrieval  uses  the 
SABER  CO2  1 5 -pm  channel  radiance  to  retrieve  Tk.  Pressure  p  is  determined 
simultaneously  with  Tk  by  using  both  narrow-filter  and  wide-filter  CO2  channels 
according  to  the  method  of  Gille  and  House  [1971],  The  retrieval  uses  two  primary 
relaxation  loops.  In  the  inner  loop,  a  Tk  profile  is  retrieved  using  the  top-down  onion-peel 
approach,  while  pressure,  CO2  vmr,  and  the  vibrational  temperatures,  Tv,b,  are  held  fixed. 
In  the  outer  relaxation  loop,  the  pressure  profile  is  rebuilt  from  bottom  up,  starting  at  the 
lower  boundary  altitude  using  the  onion-peel  retrieved  Tk  profile  and  the  barometric 
pressure  law.  The  vibrational  temperatures  are  updated  using  the  CO2  Tvib  model  with  the 
previously  retrieved  Tk,  CO2  vmr,  and  pressure  profile  as  input. 

Originally,  it  was  thought  that  CO2  was  uniformly  mixed  at  least  up  to  90  km,  so  that, 
in  the  mesosphere,  its  vmr  could  be  inserted  a  priori  into  the  forward  model  to  calculate 
temperature  [**Mertens  et  al.,  2001].  This  temperature  profile  could  then  be  inserted 
into  the  CO2  4.3-pm  forward  model  to  retrieve  CO2  vmr  from  SABER's  4.3-pm  channel. 
However,  it  became  clear  through  investigations  carried  out  under  the  SABER  program 
that  CO2  collisions  were  insufficient  to  keep  CO2  uniformly  mixed  and  as  a  consequence 
CO2  started  to  diffusively  separate  as  low  as  70  km  altitude.  This  meant  that  Tk  and  CO2 
had  to  be  retrieved  simultaneously  from  the  SABER  1 5 -pm  and  4.3-pm  channels. 

An  additional  issue  is  that  the  production  of  C02(v2)  states  radiating  at  1 5  pm  depends  on 
the  atomic-oxygen  concentration,  and  hence  so  does  the  retrieval  of  Tk  above  90  km. 
There  are  several  internal  indirect  limited  methods  of  specifying  [O]  using  SABER  data, 
but  there  is  no  routine  reliable  [O]  measurement  onboard  the  TIMED  satellite.  Hence, 
one  is  often  required  to  use  climatological  values  of  O,  despite  the  fact  that  O  is  an 
extremely  variable  species  whose  variability  is  not  inherently  predictable,  at  least  at  the 
current  time. 

The  basic  Tk  -  CO2  vmr  joint-retrieval  method  has  been  documented  in  the  literature 
[*Mertens  et  al.,  2003]  and  has  been  implemented  in  the  SABER  operational  retrieval. 
The  unique  SABER  temperature  observations  have  been  validated  [ **Mertens  et  al., 
2004]  by  comparison  with  rocketbome  falling-sphere  data  from  the  summer  MaCWAVE 
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Campaign  [ Goldberg  and  Fritts,  2004;  Goldberg  et  al.,  2004;  Fritts  et  al.,  2004]  in 
Andoya,  Norway  in  Jul  2002.  Further  validation  efforts  have  used  data  from  lidars  at 
Table  Mountain,  Calif.;  at  Ft.  Collins,  Colo.;  and  in  the  Maui/M  ALT  Campaign  at  the 
Maui  Space  Surveillance  Site  (MSSS),  also  known  as  the  AMOS  facility,  on  Haleakala, 
Maui,  Hawaii,  as  well  as  satellite  data  from  NASA’s  UARS/HALOE,  and  ESA’s 
ENVISAT-1  /  MIPAS  probes.  Examples  of  retrieved  temperature  profiles  are  compared 
to  falling- sphere  climatological  data  in  Figure  6. 


SABER-De'lvfid  Kidetic  Temperatures  far  VIA^  3  2002 
Comparison  with  Lubkon/von  Zohn  Monthly  Moon  for  March  ot  69  N 


Figure  6.  Nighttime  non-LTE  Tk  retrievals  from  SABER  measurements  on  3  March  2002  near  69N. 
The  latitude,  longitude,  orbit,  and  scan  numbers,  respectively,  for  the  SABER  measurements  are 
indicated  in  the  legend.  The  retrieved  Tk  profiles  are  compared  to  the  climatological  monthly  mean 
March  Tk  profile  from  Liibken  and  von  Zahn  [1991]  (black),  and  differences  from  the  monthly  mean 
are  shown  in  the  right  panel. 


Although  SABER  Tk  values  were  generally  in  good  agreement  with  independent 
measurements,  an  issue  arose  in  the  validation  concerning  the  location  of  the  polar 
summer  mesopause  and  the  value  of  its  associated  minimum  temperature.  In  particular,  it 
was  noticed  that  the  SABER  mesopause  was  always  up  to  3  km  lower  in  altitude  and 
colder  than  the  falling-sphere  [Liibken,  1999]  and  lidar  [Liibken  and  von  Zahn,  1991] 
high-latitude  summer  measurements.  One  consequence  was  that  noctilucent  clouds 
would  form  at  altitudes  below  80  km,  in  disagreement  with  the  results  of  ice-cloud 
microphysical  models  such  as  CARMA  ([Rapp  et  al.,  2002]  and  references  therein). 

This  issue  was  addressed  and  solved  [Kutepov  et  al.,  2006].  The  main  issue  was  that  the 
SABER  operational  forward  model  lacked  vibration-to-vibration  (V-V)  coupling  terms 
involving  redistribution  of  V2  quanta  among  CO2  isotopes.  The  AFRL  team  members 
collaborated  in  the  investigation  by  carrying  out  a  study  of  the  effect  of  V-V  coupling  on 
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CO2  V2  vibrational  temperatures  (TVI|,)  in  the  polar  summer  mesopause  region.  This  study 
involved  runs  of  AFRIT  s  ARC  (Atmospheric  Radiance  Code)  non-LTE  radiance  model 
[i Winters teiner  et  al.,  1992;  Nebel  et  al.,  1994],  This  study  verified  that  proper  treatment 
of  the  exchange  of  quanta  between  CO2  isotopes  lowers  the  summer  polar  mesopause  Tv;b 
for  minor  isotopes  by  as  much  as  25  K  while  increasing  slightly  the  major-isotope  Tv,b. 
Because  of  this,  the  discrepancies  between  SABER  retrievals  and  falling-sphere 
temperature  measurements  have  been  greatly  reduced.  Mesopause  altitude  is  in  better 
agreement,  and  retrieved  conditions  are  now  in  accord  with  noctilucent-cloud 
observations. 

Although  analysis  of  SABER  radiances  for  the  CO2  4. 3 -pm  and  15-pm  bands  have 
indicated  that  our  non-LTE  models  are  very  good  at  predicting  the  daytime  4.3-pm 
radiance,  there  is  a  considerable  underprediction  at  night.  This  is  a  longstanding  problem 
for  which  the  controversial  Kumer  mechanism  [K inner  et  al.,  1978]  was  proposed.  The 
Kumer  mechanism  involves  OH(v)  vibrational  energy  transfer  to  N2(v)  and  thence  to  CO2 
(V-V  transfer).  Since  SABER  measures  two  OH(v)  Av=l  bands  along  with  the  CO2 
bands,  it  is  capable  of  resolving  the  issue  of  whether  the  Kumer  mechanism  is  adequate 
or  a  new  mechanism  must  be  found.  Evidence  for  an  additional  nighttime  excitation 
mechanism  for  the  4.3-pm  CO2  emission  in  the  mesopause  region  has  been  published 
[**Lopez-Puertas  et  al.,  2004],  A  strong  correlation  was  found  between  C02(v3)  4.3-pm 
radiance  structures  near  85  km  and  the  strength  of  the  OH  layer  radiance  in  the  1.6-pm 
and  2.0-pm  channels.  Previous  DoD  experiments  including  ICECAP,  SPIRE,  and 
CIRRIS- 1 A  have  noted  similar  correlations  of  C02(v3)  and  OH(v)  emission.  Explaining 
the  radiance  requires  that  2.8-3  N2  vibrational  quanta  be  created  as  a  result  of  quenching 
one  OH  vibrational  quantum  by  collision  with  N2.  This  implies  a  fast,  near-gas-kinetic 
rate  for  the  V-V  transfer,  which  may  disagree  with  estimates  based  on  chemical-kinetic 
arguments. 


4.3.  Tides 

Signatures  of  the  migrating  diurnal  tide  (DT)  have  been  identified  in  SABER  CO2  15- 
pm  equatorial  radiance  data.  The  migrating  DT  is  forced  by  solar  heating.  Therefore,  it 
is  in  synchrony  with  solar  motion  and  hence  also  with  local  time.  By  separating  data 
acquisitions  into  two  halves  centered  on  the  ascending  and  descending  nodes, 
respectively,  yields  data  that  is  separated  by  approximately  nine  hours  LT  at  the  equator. 
Subtracting  this  data  gives  a  strong  DT  signal  at  the  equator  that  diminishes  gradually  as 
the  latitude  increases.  A  global  map  can  be  acquired  in  a  single  day.  By  accumulating 
such  data  for  a  full  yaw  cycle,  during  which  the  orbit  precesses  through  all  LTs 
(precession  period  ~  65  d),  the  LT  development  of  the  tidal  amplitude  and  phase  can  be 
determined  [*Picard  et  al.,  2004],  The  tidal  signature  is  present  in  both  the  CC>2(v2)  15- 
pm  radiance  and  the  retrieved  temperature  Tk.  In  Figure  7  we  see  an  example  of  the  tidal 
signal  in  Tk- 

In  the  next  two  sections,  we  discuss  the  occurrence  of  temperature-inversion  (TI) 
layers  (Section  4.4)  and  15-pm  radiance  layers  or  "knees"  (Section  4.5),  both  of  which 
turn  out  to  be  strongly  connected  with  atmospheric  tides. 
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Equatorial  Bin,  Ascending/Descending  Channel-1  Radiance  Differences 


Radiance  Difference  (per  cent  wrt  mean,  offset  2.5%  per  day) 

Figure  7.  Compilation  of  SABER  equatorial  ascending/descending  15-pm  radiance  differences, 
corresponding  to  pairs  of  LTs  separated  by  ~8.8  hours,  for  ten  days  spaced  throughout  a  33-day 
period  in  March- April  2002.  The  data  are  characterized  by  a  vertical  wavelength  of  22-25  km,  a 
growing  amplitude  with  tangent  height,  and  phase  fronts  descending  with  increasing  LT  at  a  rate 
of  ~22  km/day — all  characteristics  of  the  diurnal  tide. 


4.4.  Temperature-Inversion  Layers 


Earlier  lidar  work  (see  Meriwether  and  Gerrard  [2004]  for  a  review)  had  shown  that 
TI  layers  are  very  common  features  of  the  mesosphere  and  mesopause  region,  which  can 
be  quite  persistent  and  geographically  extensive.  The  region  above  the  TI  has  a  small 
Brunt- Vaisala  (BY)  frequency  N  and  thus  a  reduced  static  stability.  This  can  also  lead  to 


a  Richardson  number,  Ri  =  N2  / 


[du!  dxf  +(dv!  dy)2 


,  being  less  than  %,  where  u  and  v 


are  the  horizontal  wind  components  in  the  zonal  (x)  and  longitudinal  (y)  directions. 
Hence,  the  region  is  dynamically  unstable.  The  temperature  lapse  rate  above  the  TI  can 
even  be  superadiabatic,  in  which  case  the  BV  frequency  is  negative,  the  region  is 
convectively  (statically)  unstable,  and  it  will  not  support  waves. 

On  the  other  hand,  the  inversion  region  below  the  temperature  peak  has  a  large  BV 
frequency,  is  statically  very  stable,  and  can  serve  as  a  wave-duct  region.  For  further 
details  on  ducted  waves  in  TI  layers,  see  Section  5.1  below.  Earlier  work  by  the  authors 
has  shown  that  TI  layers  likely  arise  from  deposition  of  energy  and  momentum  at  a  GW 
critical  layer  [Huang  et  al.,  1998].  The  momentum  deposition  leads  to  accelerations  of 
the  background  wind,  which  can  result  in  large  wind  shears  near  the  TI  layer  and  further 
reduction  of  the  Richardson  number.  For  related  work  on  TI  layers,  see  Section  5.1.1. 
The  first  global  study  of  temperature-inversion  (TI)  layers  has  been  carried  out  under  this 
task  using  SABER  Tk  data  [** Picard  et  al.,  2004],  Inversion  layers  are  found  to  be 
globally  common  features  of  the  MLT  region,  with  amplitudes  up  to  90-100  K.  The 
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amplitudes  and  widths  of  the  layers  and  the  altitudes  at  which  the  layers  occur  exhibit 
characteristic  behavior  as  a  function  of  season,  latitude,  and  local  time  (LT).  The  latitude 
and  LT  dependence  shows  strong  DT  and  semidiurnal  tidal  (SDT)  signatures. 

An  example  of  binned  occurrence  probability  distribution  for  the  Sep-Nov  2002  yaw 
cycle  is  shown  in  Figure  8,  using  the  Level  2A  SABER  data.  The  bins  cover  2°  in 
latitude  and  lh  in  LT.  There  is  no  data  near  local  noon  since  the  sensor  cannot  look  at  the 
sun.  Considering  only  inversions  of  at  least  20  K,  we  find  TI  layers  in  80-90%  of  all 
events  in  some  regions  of  the  latitude/LT  space.  If  the  threshold  is  reduced  to  10  K,  they 
are  nearly  ubiquitous  in  these  same  regions,  which  are  mostly  nighttime  scans.  On  the 
other  hand,  in  some  other  regions,  the  chances  of  seeing  a  large  inversion  are  quite  small. 
This  is  particularly  true  during  daytime  at  high  latitude  and  midlatitude.  The  DT 
dependence  can  be  seen  near  the  equator,  while  at  northern  midlatitudes  there  is  a  hint  of 
developing  SDT  dependence.  The  tidal  dependence  is  even  clearer  in  the  binned  altitude 
distributions  (not  shown). 

The  large  spatial  extent  of  TI  layers  is  also  evident  in  the  SABER  data.  Histograms 
of  along-track  length  of  TI  layer  regions  show  that  TI  layers  are  truly  global  features, 
often  extending  for  2000-4000  km  and  more  along  the  tangent-point  track. 


TIL  Percentage,  iT>20K  Sep-Nov  2002 


Local  Time 

Figure  8.  Percentage  of  all  SABER  scans  in  which  at  least  one  inversion  layer  exceeding  20  K  was 
identified  in  retrieved  Tk,  as  a  function  of  latitude  and  LT  for  the  Sep-Nov  2002  yaw  cycle. 
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4.5.  C02 15-|nm  Radiance  Layers  /  ’’Knees” 


There  is  always  a  change  of  slope  in  the  C02  15-prn  radiance  near  110-115  km 
altitude.  This  change  has  been  shown  to  be  due  to  the  onset  of  a  new  production 
mechanism  for  the  C02(v2)  excited  state,  namely  production  by  O  collisions.  Sometimes 
the  change  of  slope  is  enough  to  produce  a  local  maximum  in  radiance.  These  intriguing 
variable  radiance  layers  or  “knees”  persistent  in  the  110-115  km  altitude  range  have  been 
characterized  on  a  global  scale  using  SABER  data.  They  are  also  found  to  exhibit  tidal 
signatures,  but  are  not  generally  associated  with  TI  layers.  Standard  radiance  models 
using  climatological  atmospheric  profiles  are  unable  to  predict  reliably  the  occurrence  of 
radiance  layers.  However,  we  have  shown  that  radiance  models  incorporating  DT 
modulation  of  the  atmospheric  profiles  are  capable  of  accounting — qualitatively  at  the 
present  time — for  the  layer  occurrence. 

Atmospheric  thermal  tides  from  the  GSWM  (Global  Scale  Wave  Model)  [Hagan  et 
al.,  1999;  http://web.hao.ucar.edu/public/research/tiso/gswm/gswm.htmll  have  been 
incorporated  into  the  ARC  non-LTE  radiance  model  and  the  results  used  to  determine  the 
influence  of  tides  on  CO2  1 5 -pm  radiance,  with  comparisons  to  SABER  data.  Although 
the  resulting  tidally  driven  model  can  account  for  the  radiance  knees,  it  is  unable  to 
account  for  the  TI  layers  which  are  equally  common  features  at  lower  altitude  in  the  MLT 
region.  By  contrast,  atmospheric  climatologies,  by  themselves,  without  inclusion  of  the 
tidal  effect  in  the  radiative  transfer,  were  unable  to  account  for  the  radiance  knees. 

4.6.  NO  5.3-nm  Radiance 

Nighttime  emissions  at  5.3  pm  from  nitric  oxide  (NO)  measured  by  SABER  were 
compared  with  the  higher-resolution  but  more  sporadic  measurements  from  the  MIPAS 
limb  emission  spectrometer  on  the  European  ENVISAT-1  probe.  The  vibrational, 
rotational,  and  spin  temperatures  were  determined  by  fitting  basis  functions  to  MIPAS 
spectral  data  [**Gardner  et  al.,  2007].  It  was  found  that  the  NO  1-0  vibrational  emission 
was  in  LTE  below  1 10  km.  Moreover,  the  rotational  and  spin  temperatures  converged 
near  110  km.  Hence,  NO  rotational  and/or  spin  temperatures  can  be  used  to  estimate 
kinetic  temperature  Tk  near  110  km.  This  is  above  the  region  where  the  standard  SABER 
temperature  retrieval  from  CO2  15-pm  emission  is  most  reliable.  These  temperatures 
near  1 10  km  were  found  to  be  considerably  above  temperatures  from  the  NRLMSISE-00 
model  [Picone  et  al.,  2002],  the  standard  semiempirical  atmospheric  model.  The  retrieval 
of  NO(v=l)  and  NO(v=0)  populations  from  the  data  were  considered.  The  latter  may 
lead  to  a  useful  retrieval  capability,  since  most  current  retrievals  (solar  fluorescence  or 
solar  occultation)  require  sunlight  and  thus  are  restricted  to  daytime  or  terminator 
conditions. 

4.7.  03  9.6-|xm  Radiance 

One  method  of  ozone  retrieval  from  SABER  data  uses  the  radiance  from  the  SABER 
9.6-pm  channel  where  O^fvj,  vi)  radiates.  However,  there  are  interfering  emissions  from 
C02(001-010)  in  this  region,  especially  in  daytime.  These  emissions  are  known  as  the 
CO2  laser  bands  and  arise  from  the  same  upper  state  that  gives  rise  to  the  first 
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fundamental  V3  band  of  CO2.  We  have  improved  the  existing  model  of  the  CO2  V3-band 
emission  by  adding  isotopic  CO2  V-V  (vibration- to-vibration)  coupling.  This  improved 
model  has  a  number  of  useful  applications.  In  particular,  it  was  used  to  resolve 
discrepancies  between  the  two  different  methods  of  retrieving  ozone  from  the  SABER 
data,  the  method  from  9.6-pm  ozone  emission  mentioned  above  and  a  method  using  the 
SABER  1 ,27-pm  CEC'Ag)  channel,  by  quantifying  the  interfering  CO2  laser-band 
emission  in  the  9.6-pm  ozone-channel  window.  This  study  showed  that  the  treatment  of 
CO2  bands  was  not  the  source  of  the  ozone-retrieval  discrepancies. 


4.8.  OH  Layers,  Interannual  Variability,  and  Tracer  Transport 

OH  Meinel-band  rovibrational  features  arise  from  a  chemiluminescent  layer  near  87 
km  altitude,  primarily  due  to  the  reaction 

H  +  03  — >  OH(v)  +  02  . 

They  serve  as  an  important  probe  of  the  mesopause  region  and  dominate  the  quiescent 
nightglow.  SABER  measures  the  Av=2  overtone  Meinel  emission  in  two  bands,  one  near 
1.6  pm  from  the  low-v  states  (v=l-3)  and  one  near  2.0  pm  arising  from  the  high-v  states 
(v=7-9). 

AFRL  6.1  team  members  have  a  long  history  of  involvement  with  OH  airglow 
modeling  (for  example,  Makhlouf  et  al.  [1995])  and  participation  in  ground-based 
measurement  campaigns,  dating  from  the  MAPSTAR,  SOAR,  and  ALOHA  programs 
[Dewan,  1994;  Makhlouf  et  al.,  1998],  TIMED/SABER  provides  the  opportunity  to 
characterize  the  variability  of  OH  Meinel  emission  on  a  global  scale  over  a  period  of 
years. 

The  SABER  OH  Meinel-band  data  have  been  used  to  determine  the  global  statistical 
characteristics  of  the  OH  layer's  altitude,  amplitude,  and  width  as  a  function  of  latitude 
and  LT.  As  an  early  finding  from  this  analysis,  OH  double,  and  even  triple,  layers  were 
found  to  be  very  common,  contradicting  the  standard  result  that  OH  was  distributed 
according  to  a  simple  Chapman-like  layer.  The  occurrence  probability  of  multiple  layers 
was  determined,  and  it  was  found  that  multiple  layers  occurred  a  very  large  fraction  (as 
much  as  40%)  of  the  time.  Double  layers  were  reported  previously  by  Melo  et  al.  [2000] 
in  analyzing  UARS/WINDII  data;  however,  the  latter  study  found  they  only  occur  5-25% 
of  the  time.  Multiple  OH-layer  occurrence  was  found  not  to  be  correlated  with  the 
existence  of  large  TI  layers. 

The  reasons  for  the  disagreement  were  investigated.  We  conjectured  that  the  seeming 
prevalence  of  multiple  layers  could  be  spurious.  Limb  images  of  narrow  layers  or 
isolated  objects  are  subject  to  inherent  ambiguity  as  to  where  along  the  LOS  structures 
giving  rise  to  bright  (or  dark)  features  are  located.  SABER  views  the  limb  perpendicular 
to  the  orbit  track,  and  uses  a  1-D  inversion  of  the  limb  radiance  profile  to  retrieve  volume 
emission  rate  (VER).  The  retrieval  uses  SABER's  isolated  slices  through  the  limb, 
resulting  in  successive  limb  profiles  about  1  minute  or  400  km  apart.  At  times  this 
produces  unrealistic  profiles  of  OH  VER,  including  peaks  at  65  km  (the  photochemistry 
of  OH  does  not  allow  for  production  of  OH  at  65  km  at  night),  double  peaks  with  a 
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secondary  peak  below  78  km,  and  negative  OH  VERs.  Artifacts  such  as  these  can  be 
produced  by  slight  variation  of  the  height  of  the  profile  along  the  LOS,  caused  by 
atmospheric  GWs. 

Work  carried  out  with  Edward  (Ted)  Llewellyn  (Univ.  Saskatchewan),  PI  for  the 
OSIRIS  instrument  [Llewellyn  et  al.,  2004]  on  the  Swedish  Odin  satellite,  showed 
[Degenstein  et  al.,  2003]  that,  as  suggested  above,  the  assumption  of  a  spherical-shell 
atmosphere  implicit  in  the  SABER  (and  WINDII)  retrievals  could  lead  to  artifacts.  The 
OSIRIS  instrument  makes  tomographic  measurements  by  scanning  in  nadir  angle  in  the 
plane  of  the  orbit  and  getting  multiple  views  through  a  radiating  volume.  Thus  it  can 
distinguish  structures  in  2-D  (altitude  and  earth-centered  angle  in  the  orbital  plane).  The 
atmosphere  was  perturbed  with  a  GW  having  a  horizontal  wavelength  of  200-400  km,  a 
large  vertical  wavelength,  and  a  vertical  displacement  of  1-2  km.  The  resulting  limb 
radiance  was  calculated  with  the  ARC  model,  and  then  the  volume  emission  rate  versus 
altitude  was  retrieved  assuming  a  spherical-shell,  horizontally  uniform  atmosphere.  This 
process  produced  double-layer  artifacts  resembling  some  of  the  double  layers  in  the 
SABER  data.  Consequently,  the  fractional  occurrence  of  double  layers  is  still  an  open 
question. 

The  SABER  OH  Meinel-band  statistics  have  been  calculated  for  the  period  2002- 
2006,  continuing  into  2007.  Interannual  variability  has  been  found  in  the  altitude  of  the 
OH  peak  at  high  northern  latitudes  in  winter.  The  high-latitude  (70-80N)  winter 
mesosphere  is  "normal"  in  2002,  2003,  and  2005,  with  an  OH  layer  peak  height  of  87-88 
km  and  an  upper  mesospheric  temperature  of  200-225  K.  On  the  other  hand,  the 
mesosphere  is  anomalous  in  2004  and  2006,  with  an  unusually  bright  OH  peak  5  km  or 
more  lower  than  the  nominal  normal  87-km  height  (81-82  km,  and  as  low  as  78  km)  and 
an  anomalously  warm  upper  mesosphere  of  260-280  K.  The  mesospheric  anomalies  are 
accompanied  by  anomalous  stratospheric  temperature  and  circulation,  with  an 
anomalously  high  stratopause  altitude  of  70-75  km  [Siskind  et  al.,  2007],  Figure  9  shows 
the  nighttime  SABER  OH  1.6-pm  band  mean  altitude  and  volume  emission  rate  statistics 
for  the  Jan-Mar  2006  yaw  period.  The  anomalously  low  bright  winter  OH  layers  are  seen 
in  the  northern-hemisphere  high  latitudes. 

These  features  can  be  explained  by  enhanced  downward  transport  of  odd  oxygen 
caused  by  the  unusual  circulation.  In  the  anomalous  periods  there  is  also  enhanced 
downward  transport  of  odd  nitrogen  (NOx)  into  the  stratosphere,  which  is  important 
because  it  results  in  destruction  of  the  ozone  layer.  Hence,  the  presence  of  very  low  and 
unusually  bright  OH  layers  in  the  mesopause  region  can  serve  as  a  warning,  or  even  a 
proxy,  for  anomalously  high  NOx  and  a  hazard  to  the  O3  layer  in  the  polar  winter 
stratosphere. 

The  analysis  of  this  data  is  continuing  into  2007.  The  year  2007  is  a  normal  year  for 
the  high-latitude  winter  mesosphere  like  2003  and  2005.  However,  the  variable  behavior 
does  not  seem  to  be  quasi-biennial  in  nature,  since  2002  was  also  a  normal  year.  Further 
correlations  with  stratospheric  dynamics  are  underway.  The  anomalous  OH  layer 
behavior  may  turn  out  to  be  a  leading  indicator  for  winter-to-spring  ozone  destruction 
events  brought  about  by  the  enhanced  downward  transport  of  NOx.  This  work  is  being 
prepared  for  publication. 


18 


Local  Time  Local  Time 


Figure  9.  SABER  nighttime  OH  1.6-jnm  band  mean-altitude  statistics  (left)  and  VER  statistics  (right)  for  the 
Jan-Mar  2006  yaw  period. 


4.9.  Gravity  Waves 

SABER  was  not  designed  to  monitor  GWs.  Yet,  using  the  technique  of  Fetzer  and 
Gille  [1994],  a  residual  method  involving  subtracting  known  zonal-mean,  planetary- 
wave,  and  tidal  perturbations  from  temperature  variance  data  and  identifying  the 
remainder  with  the  GW  variance,  temporally  smoothed  global  GW  variances  can  be 
extracted  from  a  few  consecutive  days  of  data.  The  method  was  refined  by  Peter  Preusse 
and  collaborators  [Preusse  et  al.,  2000,  2002]  who  applied  it  to  Shuttle -based  CRISTA  1 
and  2  data,  similar  techniques  were  also  applied  by  other  researchers  to  UARS  and  EOS 
Aqua  MLS  data  [Wu,  2001].  The  method  has  now  been  used  on  SABER  retrieved-Tk 
data,  giving  rise  to  global  GW  temperature-variance  maps  from  15  to  90  km  altitude. 
Using  SABER  temperature  variance  data,  the  global  characteristics  of  GWs  have  been 
determined  [** Preusse  et  al.,  2006],  The  SABER  data  represent  a  major  improvement  in 
the  altitude-range  and  temporal  coverage  using  this  method. 

Global  maps  of  GW  variance  at  various  altitudes  and  zonal-mean  latitude/altitude 
cross  sections  were  generated  from  SABER  temperature  data  for  Aug  2003  and  compared 
to  similar  data  from  the  CRISTA  instrument  on  Space  Shuttle  for  Aug  1997  and  from  the 
UARS/CLAES  instrument  for  Aug  1993  (Figure  10).  The  data  are  in  basic  agreement 
and  also  were  found  to  agree  with  runs  of  both  the  NRL  GROGRAT  GW  ray-tracer  and 
of  the  Wamer-Mclntyre  spectral  parametrization  scheme,  except  in  isolated  features 
(underestimate  of  northern  subtropical  convectively  generated  waves  and  overestimate  of 
northern  high-latitude  wave  activity  by  the  models). 

The  analysis  for  SABER  has  been  extended  to  the  first  multi-year  climatology  of 
essentially  continuous  GW  data  spanning  more  than  four  years  (2002-2006).  The 
climatology  shows  a  semiannual  variation  in  the  tropical  upper  stratosphere  and 
mesosphere  with  peaks  at  solstices,  along  with  strong  wave  activity  in  the  winter  polar 
region  and  weak  activity  near  the  summer  pole.  This  data  set  has  also  been  applied  to  the 
investigation  of  the  quasi-biennial  oscillation  (QBO)  and  semiannual  oscillation  (SAO)  in 
stratospheric  and  mesospheric  waves. 
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Figure  10.  Zonally-averaged  GW-variance  latitude-altitude  cross  sections  for  August  2003  from 
SABER  and  for  August  in  two  other  years  from  the  CRISTA-2  sensor  on  the  Space  Shuttle  and  the 
UARS/CLAES  sensor.  The  similarity  of  the  plots  stands  out,  as  well  as  the  superior  altitude  range  of  the 
SABER-derived  data. 


4.10.  Solar  Storms 

Prior  to  TIMED/SABER,  limited  measurements  of  infrared  response  to  auroral 
excitation  existed.  There  were  scattered  sounding  rocket  measurements  (for  example 
[Stair  et  al.,  1975],  a  few  measurements  from  CIRRIS-1A  [ Wise  et  al.,  1995],  and  a 
number  from  the  SPIRIT-3  sensors  on  MSX  [Sharma  et  al.,  2001;  O'Neil  et  al.,  2007], 
Although  these  earlier  experiments  provided  information  concerning  the  emitting  species 
and  helped  validate  models  of  the  non-LTE  mechanisms  they  gave  little  or  no  insight  into 
the  global  nature  of  the  response  and  the  resulting  effect  on  the  thermospheric  energy 
budget.  SABER,  with  its  virtually  1 00%  duty  cycle,  has  provided  the  first  complete  data 
and  was  fortunate  to  be  observing  during  two  major  geomagnetic  storm  periods,  April 
2002  and  October-November  2003  (the  Halloween  storm).  The  continuing 
unprecedented  solar  storms  in  the  declining  phase  of  the  current  solar  cycle  have 
provided  extraordinary  opportunities  to  document  further  the  signatures  of  emissions 
from  the  disturbed  atmosphere  in  the  IR  to  UV  spectral  range. 

Analysis  of  the  response  of  the  TIMED/SABER  satellite  IR  channels  to  the  intense 
solar-storm  period  in  April  2002,  coupled  with  observations  of  the  entire  sun-to-earth 
chain  by  approximately  thirty  other  spacecraft,  other  sensors  on  TIMED  including  the 
GUVI  uv  sensors,  and  ground-based  instruments,  has  yielded  the  most  detailed 
knowledge  to  date  of  our  understanding  of  the  upper-atmospheric  radiance  response  to 
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solar  storms.  AFRL  6.1  team  personnel  were  involved  in  analyzing  and  understanding 
the  wavelength  channel  dependence  of  the  storm  response.  Detailed  descriptions  of 
energy  transport  in  the  Apr  2002  solar  storm  based  on  SABER  data  have  been  published 
\**Mlynczak  et  al.,  2003;  Winick  et  al. ,  *2004a,  *2004b;  **Mlynczak  et  al.,  2005], 
documenting  the  flow  of  energy  into,  out  of,  and  within  the  thermosphere.  To  study  the 
energy  flow,  the  SABER  data  was  enhanced  by  runs  of  the  ASPEN  thermospheric 
general  circulation  model.  Enhancements  of  factors  of  thirty  were  observed  in  the 
SABER  NO  channel  at  5.3  pm,  accompanied  by  rich  structure  with  strong  latitudinal  and 
longitudinal  gradients.  Auroral  energy  was  transported  into  low  latitudes  and  the 
equatorial  region  over  a  time  scale  of  two  to  three  days  and  was  accompanied  by  sharp 
increases  in  low  and  mid-  latitude  radiance.  This  was  followed  by  reestablishment  of 
quiescent  conditions  over  the  course  of  a  few  days  by  a  dramatic  increase  in  the  NO  5.3- 
pm  cooling  rate  to  -2200  K/d,  a  negative-feedback,  stabilizing  process  which  we  have 
called  the  “nitric  oxide  thermostat”.  Reradiation  by  NO  dissipated  50%  of  the  energy 
input  to  the  thermosphere  by  the  Apr  2002  storm. 

The  team  has  also  investigated  storm  effects  in  various  other  spectral  channels.  These 
include  the  CO2  4.3-pm  band,  the  O2C 1  Ag)  1.27-pm  band,  the  OH  2.0-pm  band,  and,  to  a 
lesser  extent,  the  OH  1 ,6-pm  bands.  Some  of  these  bands,  in  particular  4.3-pm,  show 
order-of-magnitude  enhancements.  The  enhancement  in  the  4.3-pm  channel  is  not 
entirely  due  to  CO2,  but  includes  emission  from  aurorally  enhanced  NO+  ions,  discussed 
further  in  the  next  section.  The  SABER  4.3-pm  channel  radiance  enhancement  maps 
show  good  correlation  with  similar  maps  of  far-uv  LBH  (Lyman-Birge-Hopfield)  band 
emission,  clearly  auroral,  from  the  GUVI  sensor  onboard  TIMED.  Note,  however,  also 
that  the  solar-storm  radiators  in  the  OH  channels  are  not  related  to  OH.  There  has  been 
some  inconclusive  investigation  of  the  nature  of  these  radiators;  one  strong  candidate  for 
the  2.0-pm  channel  is  NO+-ion  overtone  emission. 

In  addition,  the  team  investigated  the  response  of  SABER  emissions  to  solar  flares 
[*  Winick  et  al.,  2004b],  The  flares  are  short-duration  daytime  events  which  affect  all 
latitudes,  even  equatorial.  Because  of  the  TIMED  orbit,  it  may  take  up  to  12  hours  for 
the  SABER  instrument  to  view  the  affected  region.  To  minimize  solar- wind 
(geomagnetic)  effects,  it  is  best  to  look  at  low  magnetic  latitude  after  a  period  of  low 
geomagnetic  activity  (low  Kp).  The  SABER  data  showed  that  there  was  a  significant 
response  of  the  NO  5.3-pm  radiance  to  the  historically  large  X18  solar  X-ray  flare  of  28 
Oct  2003,  just  before  the  intense  Halloween  solar  storm.  Figure  1 1  shows  the  clear 
response  to  the  28  Oct  flare  in  the  SABER  5.3-pm  band  radiance. 


4.11.  NO+(v)  4.3-pm  Emission,  Thermospheric  Temperature,  and  E- 
Region  Solar  Storm  Response 

After  identifying  NO+  in  auroral  spectra  from  the  rocketbome  FWI  experiment 
[, Picard  et  al.,  1987],  the  AFRL  team  showed  quantitatively,  based  on  CIRRIS-1A  data, 
that  rovibrational  emission  from  the  NO+  ion  is  a  significant  contributor  to  radiance  in  the 
4.3-pm  channel  under  daytime  or  auroral  conditions,  in  addition  to  the  main  CO2  V3 
emission.  This  is  significant  for  SABER  since,  as  mentioned  above,  this  channel  is  used 
to  derive  temperature  and  CO2  vmr.  Hence,  solar  storms  or  photoelectrons  generated  by 


21 


October  28  Flare  Effect  on  NO  5,JLm  Radiance  Low  Latitude 


solar  uv  can  enhance 
the  emission  in  this 
channel  and  produce 
biases  in  the  SABER 
temperatures,  since  Tk 
is  retrieved  in  the 
SABER  operational 
software  ignoring  NO+ 
emission  in  the  4. 3 -pm 
channel.  On  the 
positive  side,  there  is  a 
potential  for  global 
measurements  of  4.3- 
pm  emission  to  give 
insight  into  E-region 
ion-neutral  chemistry 
and  energy  transfer. 

The  response  of  the  E- 
region  to  auroral 
electron  dosing  was 
investigated,  based  on 
global  measured 
SABER  nighttime  4.3- 

pm  emission  and  SABER  retrieved  temperature,  pressure,  and  CO2  vmr  [**Mertens  et 
al.,  2007].  This  study  provides  a  unique  perspective  on  E-region  chemistry,  energetics, 
radiative  transfer,  and  response  to  solar  geomagnetic  storms.  By  modeling  and 
subtracting  the  emission  due  to  CO2  from  the  measured  4.3-pm  radiance,  one  is  left  with 
NO+  emission,  which  can  be  related  to  E-region  electron  dosing  through  the  FLIP  model 
[ Richards ,  2002],  The  SABER  NO+  VER  was  shown  to  be  a  useful  proxy  for  the  storm¬ 
time  E-region  response  and  will  be  used  to  develop  parametrizations  for  storm-time 
corrections  to  the  IRI  (International  Reference  Ionosphere)  [ Bilitza ,  2001]  electron  and 
ion  densities.  Analyses  have  been  carried  out  for  the  April  2002  and  Halloween  2003 
storm  periods.  6. 1  team  member  Jeremy  Winick  and  NASA  Langley  investigator  Chris 
Mertens  are  continuing  this  work  with  partial  funding  from  a  NASA  Office  of  Space 
Science  Guest  Investigator  Grant. 

For  discussion  of  NO+  emission  in  the  MSX  data,  see  Section  3.3  above. 


Figure  11.  NO  5.3-pm  emission  observed  by  SABER  before  the 
XI 8  solar  X-ray  flare  of  28  Oct  2003  and  after  the  flare  at  low 
latitude.  Each  curve  shows  the  zonal  mean  for  the  day,  while  the 
bars  show  the  extent  of  the  diurnal  variation. 


4.12.  Auroras 

Possible  auroral  monitors  internal  to  SABER  data  at  1.6  and  2.0  pm  have  been 
identified.  The  2.0-pm  auroral  contribution  is  the  stronger  of  the  two.  SABER  was 
explicitly  designed  to  avoid  auroral  emissions.  However,  these  channels  that  were 
intended  to  monitor  chemiluminescent  OH  Meinel-band  emission  at  tangent  heights  of 
83-96  km  show  a  response  at  high-tangent  heights  >100  km  that  is  correlated  with  auroral 
energy  deposition.  This  is  particularly  opportune  since  other  TIMED  instruments  (GUVI 
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and  TIDI)  that  could  provide  auroral  inputs  to  SABER  have  proven  to  be  less  than  ideal 
for  this  purpose,  resulting  in  team  members  devoting  considerable  time  to  looking  for 
internal  auroral  monitors.  Work  is  continuing  on  identifying  the  source  of  the  emission 
and  its  properties.  Some  candidates  are  NO+  and  possibly  N2+  ionic  emissions  [*D.  R. 
Smith,  2003],  The  non-auroral  background  in  these  channels  above  100  km  tangent 
height  is  very  small,  since  one  is  above  the  OH  emission  layer. 

Other  SABER  spectral  bands  also  respond  to  auroral  dosing  [*  Winick  et  al.,  2004a], 
as  discussed  in  Sections  4.10  and  4.1 1.  (1)  The  NO  5.3-pm  emission  from  the  SABER 
5.3-pm  channel  has  an  airglow  component  in  the  1-0  band  that  is  enhanced  during  aurora. 
In  addition,  there  is  copious  production  of  hot  fundamental-band  emissions  from  NO 
during  auroras  [ Winick  et  al.,  1987].  However,  the  fundamental  auroral  hot  bands  fall 
largely  outside  the  SABER  5.3-pm  channel  bandpass.  As  discussed  in  Section  3.3, 
vibrationally  excited  NO  is  produced  immediately  by  chemiluminescent  reactions 
between  N(2D)  and  O2.  This  process  will  also  enhance  total  NO.  The  NO  produced 
could  then  be  excited  by  atomic  oxygen  (the  prime  high-altitude  quiescent  NO(v=l) 
production  mechanism),  leading  to  enhanced  airglow  and  enhanced  radiance  in  the 
SABER  5.3-pm  channel.  (2)  The  C02(v3)  4.3-pm  emission  has  a  strong  nighttime 
quiescent  contribution,  but  is  enhanced  in  a  strong  aurora  (IBC  III  or  greater)  [Stair  et  al., 
1975;  Winick  et  al.,  1987].  The  NO+  rovibrational  fundamental  also  radiates  near  4.3  pm, 
and  at  night  it  is  largely  of  auroral  origin.  Both  the  CO2  and  NO+  emission  will  enhance 
the  radiance  in  the  SABER  4.3-pm  bandpass  during  auroras.  (3)  Since  the  CO2  laser 
bands  at  9.4  and  10.4  pm  also  originate  from  the  same  V3  upper  state  as  the  4.3-pm 
emission,  they  are  also  enhanced.  This  will  increase  the  amount  of  interference  from  CO2 
in  the  SABER  O3  9.6-pm  channel.  (4)  The  Cfl'Ag)  1.27-pm  band,  which  is  measured  by 
the  shortest-wavelength  SABER  radiometer  channel,  is  enhanced  during  aurora.  This 
emission  is  most  likely  due  to  excitation  of  O2  by  auroral  electrons. 

We  examined  the  correlation  between  the  SABER  channel  and  TIMED/GUVI  LBH-2 
channel  radiance — a  good  proxy  for  the  auroral  dosing  [*  Winick  et  al.,  2004a],  It  was 
found  that  the  correlation  was  the  strongest  for  the  2.0-pm  channel,  which  appears  to 
have  a  strongly  linear  response  to  dosing  and  the  purest  auroral  response.  On  the  other 
hand,  the  NO  5.3-pm  channel  correlation  is  the  weakest,  since  there  is  a  large  quiescent 
component.  In  addition,  the  auroral  response  may  be  nonlinear  and  nonlocal,  since  NO  is 
produced  by  auroras,  has  a  long  lifetime,  and  is  advected  by  winds.  The  correlation  of 
the  SABER  4.3-pm  channel  with  the  GUVI  LBH  radiance  is  intermediate  between  the 
two  above.  Once  again,  there  is  a  significant  quiescent  component.  However,  at  higher 
dosing  levels  the  airglow  component  is  less  important,  there  is  less  scatter,  and  a  better 
linear  correlation  is  established. 


5.  RADIANCE  STRUCTURE 

The  nonequilibrium  middle  and  upper  atmosphere  (30-160  km)  are  very  dynamic 
regions  that  are  structured  vertically  and  horizontally  by  numerous  processes.  These 
include  GWs,  turbulence,  fronts,  layers,  auroras,  the  terminator,  thundercloud  E  fields, 
upward  lightning  (sprites  and  their  kin),  thermal  plumes,  and  meteor  trails.  Of  particular 


23 


interest  is  the  passage  of  both  atmospheric  GWs  and  transient  frontal  disturbances  or 
bores,  as  well  as  the  presence  of  persistent  TI  layers.  The  infrared  emissions  from  this 
part  of  the  atmosphere,  sometimes  called  the  "ignorosphere",  are  not  well  characterized. 
They  are  typically  not  in  LTE  with  the  atmosphere  and  are  further  perturbed  by  the 
presence  of  the  mentioned  sources  of  pervasive  atmospheric  structure.  The  inevitable 
result  is  highly  structured  atmospheric  emissions  that  depend  on,  but  do  not  simply 
mirror,  the  structure  of  the  atmosphere. 

Understanding  the  structure  of  the  atmosphere  is  essential  to  understanding  the 
structure  of  the  radiation  that  it  emits.  At  the  same  time,  understanding  how  atmospheric 
structure  perturbs  atmospheric  radiation  provides  a  means  to  sense  the  perturbing 
atmospheric  processes  remotely. 

We  have  examined  methods  to  calculate  the  LTE  /  non-LTE  radiative  response  to 
temporal  and  spatial  variations  of  the  atmosphere  and  given  examples  of  applications. 
Furthermore,  we  have  compared  our  results  with  existing  field  data.  In  addition,  we  have 
found  new  phenomena  in  the  underlying  atmosphere  (for  example,  undular  bores),  and 
more  completely  described  known,  but  inadequately  characterized,  phenomena  (for 
example,  GW  power  spectra).  We  have  studied  the  features  and  properties  of  these 
phenomena  independently  from  the  optical/IR  radiation  they  emit,  and,  in  many  cases,  we 
have  also  examined  the  resulting  radiative  structure.  Further,  the  team  has  participated  in 
the  formulation  of  a  proposed  new  and  as  yet  unfunded  NASA  optical  /  infrared 
experiment  (Waves  Explorer),  whose  goal  is  to  sense  atmospheric  GWs  remotely  from 
earth  orbit  on  a  global  basis  and  characterize  their  sources.  We  discuss  localized  sources 
of  structure  in  Section  5.1  and  structuring  by  GWs  and  other  phenomena  in  Section  5.2. 


5.1.  Localized  structures 

Atmospheric  structures  and  the  radiance  structures  they  produce  can  be  localized  in 
the  vertical  direction,  or  in  the  horizontal  direction,  or  both.  TI  layers  are  vertically 
localized  structures  having  a  strong  effect  on  emitted  and  scattered  radiation.  These  TI 
layers  will  typically  be  accompanied  by  radiance  layers.  In  addition,  in  Sections  3.2  and 
4.4  above,  we  have  described  radiance  layers,  or  knees,  that  have  a  more  complex  origin 
and  are  not  associated  with  TI  layers.  Temperature  structure  in  general  and  TI  layers  in 
particular  can  produce  GW-ducting  regions.  These  ducting  regions  confine  GWs  and 
localize  them  in  the  vertical  direction.  Such  wave  ducts  can  also  be  formed  by 
horizontal- wind  curvature  (Doppler  ducting)  [ Chimonas  and  Hines,  1986]  or  by  a 
combination  of  both  temperature  and  wind  structure  [Fritts  and  Yuan,  1989], 

Horizontally  localized  structures  of  importance  include  plumes,  jets,  and  especially 
fronts.  Fronts  can  involve  primarily  mass  flow  (gravity  currents)  or  primarily  phase 
motion  (wave  fronts)  [Simpson,  1997].  In  earlier  work  we  discovered  a  new  source  of 
frontal  disturbances  in  the  MLT  called  internal-wave  mesospheric  undular  bores.  These 
wavefronts  propagate  in  GW  ducts  and  are  seen  as  dramatic  features  of  airglow  imagery. 
Vertical  localization  due  to  layered  radiating  minor  species  is  also  possible.  Such  layers 
have  been  seen  in  the  stratosphere  in  the  form  of  narrow  ozone  laminae  (for  example, 
[Reid  and  Vaughan,  1993]).  Other  minor-species  layers  have  been  seen  in  tropospheric 
and  stratospheric  differential-absorption  lidar  profiles.  This  type  of  layering  usually 
involves  dynamical  processes  along  with  large  gradients  of  minor-species  vmr.  The 
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MLT  does  have  several  common  layered  minor  species,  occurring  in  Chapman-like 
layers,  whose  layering  is  due  to  chemistry.  Photochemistry  results  in  layers  that  are  quite 
broad,  this  layering  is  generally  well  understood.  Examples  are  the  Na  and  OH  layers. 
However,  to  our  knowledge,  narrow  minor-species  layers  of  dynamical  origin  have  not 
been  reported  in  the  MLT  region,  and  we  will  not  discuss  them  further.  However,  such 
layers,  if  present,  would  be  expected  to  lead  to  enhanced  radiance  layers  in  a  similar  way 
to  TI  layers. 

In  Section  5.1.1  we  discuss  TI  layers  further.  Section  5.1.2  deals  with  further 
comments  on  radiance  layers,  and  Section  5.1.3  discusses  mesospheric  fronts,  bores,  and 
wave  ducts. 

5.1.1.  Temperature-Inversion  Layers 

A  case  study  on  a  mesospheric  TI  layer  seen  by  the  Univ.  of  Illinois  sodium  wind- 
temperature  lidar  in  the  NSF/AFRL  ALOHA  93  Campaign  at  the  AF  MSSS  site  (AMOS) 
has  been  revisited  [**Huang  et  ah,  2002],  showing  conclusively  that  the  TI  layer  resulted 
from  GW  critical-layer  interaction.  GW  critical-layer  interactions  occur  when  the 
horizontal  phase  speed  of  the  wave  matches  the  background  wind  speed.  When  a  critical 
layer  occurs,  the  wave  can  interact  with  the  background  flow  continuously  and  exchange 
energy  with  it.  This  can  result  in  extraction  of  energy  from  the  wave  and  acceleration  of 
the  background  wind.  The  wind  acceleration  can  lead  to  a  region  of  large  shear  which  is 
dynamically  unstable  (Ri  <  0.25).  In  the  resulting  region  of  dynamic  or  shear  instability, 
the  wind  shear  produces  turbulence  and  eventually  heat,  raising  the  temperature  in  a 
narrow  TI  layer. 

The  nature  and  source  of  these  prevalent  and  persistent  layers  that  provide  stressing 
backgrounds  for  IR  surveillance  systems  has  been  controversial.  The  original  publication 
on  the  critical-layer  mechanism  for  TI  layers  [Huang  et  al.,  1998],  also  involving  team 
members,  made  the  case  for  it,  but  used  a  less  complete  version  of  the  ALOHA  lidar  data 
set  in  calculating  the  BV  frequency.  This  showed  the  existence  of  dynamically  unstable 
regions  (Ri  <  %).  The  reanalysis  shows  that  there  are  also  numerous  regions  of 
convective  instability  in  the  data,  where  N2  is  negative  and  Ri  <  0.  These  regions  also 
lead  directly  to  turbulent  heating.  Thus  both  shear  and  convective  instability  play 
important  roles  in  the  TI  layer  development.  The  reanalysis  explains  new  aspects  of  the 
observations,  including  the  lack  of  local  observation  of  the  implicated  GW  in  the  MSSS 
CCD  airglow  imagery. 

Satellite  databases  provide  ideal  opportunities  to  study  TI  layers  on  a  global  basis.  A 
global  analysis  of  the  SABER  data  for  TI  layers  is  described  in  Section  4.4. 


5.1.2.  Radiance  Layers  /  Knees 

We  have  mentioned  above  that  radiance  layers,  or  knees,  seem  to  have  a  complex 
origin.  The  striking  layers  in  1 5 -pm  CO2  emission  between  100  and  115  km  can  occur 
without  the  presence  of  any  obvious  TI  layers.  They  appear  to  involve  an  elaborate 
interaction  between  atmospheric  tides  and  radiative  excitation/emission  processes. 

Unlike  TI  layers,  such  radiance  structures  are  difficult  or  impossible  to  observe  from 
the  ground  due  to  a  combination  of  absorption  by  the  lower  atmosphere  and  foreground 
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path  radiance  resulting  from  lower-atmospheric  emission.  However,  satellite  remote¬ 
sensing  databases,  such  as  offered  by  MSX  and  SABER,  are  ideal  tools  for  investigating 
the  global  behavior  of  radiance  layers.  The  findings  from  MSX  and  SABER  on  radiance 
layers  were  discussed  in  Sections  3.2  and  4.5. 

5.1.3.  Fronts  /  Mesospheric  Bores  and  Ducted  Waves 

A  striking  front  covering  the  entire  sky  and  followed  by  phase-locked  waves  was 
observed  by  ground-based  CCD  mesospheric-airglow  imagers  [ Taylor  et  al.,  1995]  taking 
part  in  the  ALOHA  93  Campaign  from  the  AMOS  site  on  Haleakala.  The  front  was 
dubbed  the  "Great  Wall".  When  the  sharp  front  was  crossed,  dramatic  changes  in 
radiance,  either  positive  or  negative,  occurred.  Such  steps  in  radiance  would  result  in 
stressing  backgrounds  for  surveillance  systems.  AFRL  scientists  had  explained  the  front 
[Dewan  and  Picard,  1998]  as  resulting  from  propagation  of  an  undular  bore  in  an 
intemal-GW  duct.  Such  bores  are  nonlinear  propagating  phase  fronts  and  are  well- 
known  in  many  river  estuary  channels,  as  in  the  famous  Bay  of  Fundy  bore 
rhttp://www. centralnovascotia.com/tides.phpT  They  have  also  been  observed  in  oceanic 
internal  waves  [ Osborne  and  Burch,  1980]  and  in  waves  on  the  atmospheric  boundary 
layer.  The  latter  include  the  so-called  "morning  glory"  in  the  Gulf  of  Carpinteria, 
Australia  [Smith,  1988]  and  the  bores  seen  in  the  Central  Plains  in  thunderstorm  outflows 
[Doviak  et  al.,  1989],  Dewan  and  Picard  [1998]  adapted  the  existing  theory  for  channel 
bores  [Lighthill,  1978]  to  internal-wave  bores  by  requiring  that  the  internal  bores  possess 
a  nodal  line  in  the  center  of  the  duct.  These  so-called  varicose  modes  result  in 
geometrically  symmetric  patterns,  while  the  resulting  vertical  displacements  are 
described  by  antisymmetric  functions.  The  nodal  line  plays  the  role  of  the  rigid  ground 
boundary  in  the  channel  bores.  The  sharp  wave-structured  fronts  associated  with 
mesospheric  bores  should  lead  to  large  gradients  of  background  radiance  seen  by  IR 
surveillance  systems. 

Subsequent  to  this  work,  a  small  explosion  of  activity  in  airglow  studies  of 
mesospheric  bore  and  frontal  propagation  occurred.  It  was  found  that  bores  were 
infrequent,  but  not  rare,  events  seen  by  many  observers  in  all  latitudes  and  seasons. 
Researchers  have  studied  their  occurrence  characteristics  at  various  geographically 
diverse  sites.  For  examples,  see  [*Taylor  et  al.,  2003],  Surprisingly,  according  to  some 
researchers,  bores  may  occur  as  much  as  10%  of  the  time  in  the  equatorial  zone  [Fechine 
et  al.,  2005], 

Natural  questions  that  arise  are  what  causes  mesospheric  bores  and  how  are  they 
generated.  The  first  paper  on  the  generation  mechanism  for  mesospheric  undular  bores, 
was  published  by  ** Dewan  and  Picard  [2001].  It  was  suggested  that  the  body  force 
required  to  generate  the  bore  was  provided  by  a  GW  exchanging  momentum  and  energy 
with  the  atmosphere  in  a  critical-layer  interaction.  As  discussed  in  Section  5.1.1,  GW 
critical-layer  interactions  were  proposed  [Huang  et  al.,  1998,  **2002]  as  a  mechanism 
for  producing  TI  layers,  which  can  provide  regions  of  increased  stability  capable  of 
ducting  waves.  **Dewan  and  Picard  [2001]  proposed  that  it  is  these  very  same  critical- 
layer  interactions  that  can  launch  bores.  The  body  force  generated  at  the  critical-layer 
acts  in  piston-like  fashion  to  produce  a  perturbation,  which  eventually  becomes  a  bore  by 
nonlinear  steepening  as  the  perturbation  propagates  through  some  distance.  The  method 
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of  characteristics  is  used  in  this  paper  to  predict  the  time  required  to  form  the  bore,  which 
is  in  agreement  with  data. 

In  an  appendix  this  paper  also  enumerates  the  properties  by  which  an  airglow  frontal 
observation  can  be  identified  as  a  mesospheric  bore.  The  chief  characteristic  of  note  is 
the  intensity  step  at  the  front.  Such  a  sharp  intensity  step  cannot  arise  out  of  linear  GW 
interactions,  since  rectifying  a  non-zero-frequency  GW  requires  that  one  be  in  the 
nonlinear  wave  propagation  regime. 

The  first  simultaneous  and  collocated  observations  of  a  mesopause-region  undular 
bore  and  of  the  conjectured  TI  layer  providing  a  duct  for  it  were  reported  [**She  et  al., 
2004],  The  Colorado  State  Univ.  (CSU)  sodium  wind-temperature  lidar  at  Ft.  Collins, 
Colorado  made  observations  of  the  TI  layer  and  the  wind  shift  associated  with  the  bore, 
while  the  CSU  OH  airglow  CCD  imager  at  Platteville  observed  the  horizontal  structure  of 
the  bore.  Comparisons  were  made  with  the  Dewan-Picard  simple  analytic  bore  model 
[Dewan  and  Picard,  1998,  2001],  extended  to  include  background  wind.  In  this  same 
publication  the  transition  of  an  undular  bore  to  a  turbulent  (or  “foaming”)  bore  was 
observed  for  the  first  time  (Figure  12).  The  bore  must  dissipate  energy,  but  the  foaming 
bore  dissipates  energy  through  its  turbulent  front  rather  than  by  shedding  a  short-period 
wave  train.  The  foaming  bore  is  more  difficult  to  identify  through  airglow  data  than  the 
undular  bore  because  many  observers  routinely  produce  difference  images  (between 
successive  frames)  from  their  CCD  airglow  imagery,  in  order  to  enhance  signal-to-noise 
ratio  for  wave  perturbations.  This  means  that  the  uniform  background  level  change 
across  the  bore  front  will  not  be  seen  and  the  only  identifying  characteristic  of  the  bore 
remaining  will  be  a  localized  perturbation  at  the  frontal  boundary.  The  Richardson 
number  Ri  profile  is  calculated  from  the  CSU  lidar  wind  and  temperature  data,  and  it 
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Figure  12.  OH  Meinel  CCD  difference  images  (left)  showing  development  of  undular  bore 
into  turbulent,  or  foaming,  bore  over  Ft.  Collins,  Colorado  on  the  night  of  6/7  Oct  2002  near 
07:00  UT.  On  the  right  the  Richardson  number  Ri  calculated  from  the  CSU  lidar  wind  and 
temperature  data  shows  that  Ri  suddenly  drops  below  0.25  at  OH  emission  altitudes 
simultaneously  with  the  bore  becoming  turbulent. 
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shows 

(Figure  12)  that  the  atmosphere  suddenly  becomes  dynamically  unstable  (Ri  <  %)  in 
the  region  of  OH  emission  at  the  very  time  the  undular-to-turbulent  transition  occurs. 

In  addition  to  the  sharp  front,  the  amplitude  change  across  the  front,  the  characteristic 
propagation  speed,  and  the  associated  waves,  the  sense  of  the  amplitude  change  across 
the  front  is  significant.  Taylor  et  al.  [1995]  were  recording  four  airglow  emission  layers 
at  the  time  of  the  ALOHA  93  bore.  They  found  that  the  lower-altitude  OH-Meinel  and 
Na  D-line  emissions  brightened  as  one  crossed  the  front,  corresponding  to  downward 
displacement  of  the  air  parcels,  while  the  higher-altitude  0( 1 S)  (or  01)  green- line  and 
02(1Z)  Atmospheric-band  emissions  became  less  bright,  corresponding  to  upward 
displacement.  In  the  varicose-mode  bore  interpretation  of  Dewan  and  Picard  [1998],  this 
situation  corresponds  to  a  nodal  line  occurring  between  the  Na  and  O2  Atmospheric 
layers.  We  will  call  this  a  "normal"  bore,  since  the  nonlinear  steepening  that  forms  the 
bore  can  only  occur  if  the  displacement  of  the  parcels  as  the  bore  passes  is  positive.  This 
follows  from  the  fact  that  the  phase  speed  increases  with  amplitude.  Negative  bores, 
otherwise  known  as  "inverse"  bores,  are  not  stable  and  will  tend  to  flatten  and  disappear 
as  they  propagate. 

A  large  database  of  two-color  observations  from  Utah  State  University's  Bear  Lake 
Observatory  has  been  searched  for  bores  and  related  frontal  phenomena  [*  Taylor  et  al., 
2003].  These  bores  have  been  described,  analyzed,  and  compared  to  bore  theory.  Many 
examples  of  normal  bores  were  found.  Some  of  these  had  complementary  behavior  of 
the  low-altitude  (OH)  emission  and  the  high-altitude  (01)  emission,  with  a  contrast 
reversal  between  the  emissions,  just  like  the  ALOHA-93  bore.  Others  had  no  contrast 
reversal,  with  both  emissions  increasing  in  brightness  across  the  front.  In  the  varicose 
bore  model,  this  corresponds  to  a  nodal  line  occurring  above  both  of  the  emitting  layers. 

It  was  not  long  before  observers  found  examples  of  "impossible"  inverse  bores.  The 
CARI  (Cornell  All-sky  Roving  Imager)  system  found  inverse  bores  in  measurements  at 
Arecibo,  Puerto  Rico  on  6  Mar  1997  and  at  Colan-Piura,  Peru  on  24  Dec  2000.  In  these 
examples  a  bright  front  occurred  in  the  high-altitude  emission  above  a  dark  front  in  the 
low-altitude  emission — apparently  caused  by  a  negative  bore  that  is  impossible  according 
to  bore  theory.  We  suggested  that  this  did  not  represent  an  observation  of  an  inverse  bore 
but  rather  indicated  a  dual  bore,  with  one  bore  nodal  plane  above  the  higher  airglow  layer 
and  the  second  below  the  lower  layer  [*Loughmiller  et  al.,  2006]  .  The  two  bore 
responses  are  synchronized,  but  there  is  a  horizontal  phase  lag  in  the  upper-layer  bore 
data  consistent  with  the  propagation  time  for  an  upward-propagating  GW  of  appropriate 
wavelength  between  the  two  layers.  This  suggests  that  a  common  upward  GW  generated 
ducts  for  both  waves.  The  ducts  resulted  from  two  TI  layers  formed  by  critical-layer 
interaction  with  a  background  wind  modulated  by  the  diurnal  tide.  The  vertical 
wavelength  of  the  tide  is  ~20  km  which  is  greater  than  the  separation  of  the  two  airglow 
layers  in  question,  so  the  critical  layers  involved  could  have  been  formed  by  GW-tidal 
interactions  at  successive  antinodes  of  the  tidal  horizontal  velocity  [ *Loughmiller  et  al., 
2006], 

Observations  of  OH,  sodium  and  atomic-oxygen  green-line  airglow  at  the  AMOS 
site  by  the  CASI  (Cornell  All-Sky  Imager)  instrument  as  part  of  the  AFOSR-NSF  Maui- 
MALT  program  were  analyzed  for  the  occurrence  of  bores  [*Loughmiller  et  al.,  2007]. 
Nine  candidates  were  found  in  the  period  Jun  2002  -  Jun  2003,  including  examples  of 
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both  normal  bores  and  dual  (inverse)  bores.  The  imager  data  was  correlated  with  SABER 
space-based  temperature  structure  data  from  near-coincident  limb  scans.  A  number  of  TI 
layers  were  found  in  the  SABER  data,  resulting  in  the  presence  of  Brunt  (thermal)  wave 
ducts  and  corroborating  the  bore  types  seen  in  the  imager. 

On  the  theoretical  front,  progress  has  been  made  in  developing  a  theory  and 
numerical  model  for  mesospheric  undular  bores  that  allows  exploring  the  relations 
between  such  bores,  nonlinear  wave  interactions,  and  solitons.  The  theory  proposed  by 
Dewan  and  Picard  [1998],  while  extremely  useful,  contains  numerous  simplifying 
assumptions  about  the  atmosphere  and  does  not  apply  to  general  atmospheres  nor  permit 
detailed  calculations.  Moreover,  it  is  not  completely  internally  consistent  [ Lighthill , 

1978,  pg  466-467]  and  does  not  make  connections  with  the  larger  body  of  nonlinear  wave 
theory.  On  the  other  hand,  the  numerical  model  based  on  nonlinear  wave  theory  (NWT) 
allows  better  understanding  of  the  requirements  for  the  existence  of  bores  under  more 
general  conditions  and  better  prediction  of  their  occurrence.  The  NWT  approach  is  less 
general  than  direct  numerical  simulation  (DNS);  however,  its  results  are  more  easily 
interpreted  than  DNS  results,  since  DNS  is  basically  a  numerical  experiment  and 
interpretation  of  its  results  requires  the  same  degree  of  ingenuity  as  interpreting  the 
results  of  an  experiment.  The  NWT  model  was  developed  for  the  asymptotic  long-wave 
regime  and  leads  to  clear  interpretations  within  its  range  of  validity.  In  the  NWT  the 
solutions  factor  and  give  rise  to  two  coupled  equations.  These  are  the  Taylor-Goldstein 
(TG)  equation  for  the  vertical  behavior  of  the  field  and  the  Benjamin-Davis-Ono  (BDO) 
nonlinear  wave  equation  for  the  horizontal  and  temporal  behavior  in  the  deep-fluid  case. 
The  shallow- fluid  case  leads  to  a  different  nonlinear  wave  equation,  the  Korteweg- 
deVries  (KdV)  equation 

Results  from  the  NWT  are  currently  being  prepared  for  publication.  These 
include  the  following:  (1)  Bores  require  a  stable  region  in  which  to  propagate. 

Unducted  free  waves  cannot  form  bores.  (2)  The  stable  region  required  by  bores  does  not 
necessarily  require  an  inversion  layer.  Regions  of  wind  curvature  can  also  contribute  to 
duct  formation.  NWT  predicts  bores  can  exist  in  Doppler  ducts  as  well  as  in  Brunt  ducts 
formed  by  inversion  layers.  The  basic  requirement  is  that  there  be  a  region  where  the 
Scorer  parameter  is  positive,  bounded  by  regions  of  zero  or  negative  Scorer  parameter. 

(3)  Only  positive  bores  are  stable.  (4)  The  undulations  behind  the  wavefront  grow  in 
number  and  increase  in  amplitude  as  the  bore  propagates.  (5)  The  NWT  shows  clearly 
that  the  stability  of  bores  results  from  balancing  dispersion  with  nonlinear  steepening. 
While  infinitesimal-amplitude  waves  in  a  shallow  fluid  are  dispersionless  to  lowest  order, 
waves  in  a  narrow  waveguide  (that  is,  a  waveguide  whose  thickness  is  much  less  than  the 
horizontal  wavelength)  embedded  in  a  deep  fluid  are  dispersive.  (6)  Small-amplitude 
bores  travel  at  the  linear  long-wave  phase  speed  Co,  while  large-amplitude  bores  travel  at 
speeds  greater  than  co.  (6)  The  waveforms  of  small-amplitude  bores  look  like  damped 
sinusoids,  whereas  those  of  large-amplitude  bores  are  pulse-train-like  and  look  distinctly 
cnoidal,  with  sharp  crests  and  broad  flat  troughs.  Individual  pulses  are  starting  to  look 
like  solitons.  (7)  The  NWT  equations  support  soliton  solutions,  where  the  balance  of 
nonlinearity  and  dispersion  is  perfect.  For  the  KdV  equation,  the  solitons  take  the 
familiar  sech2  form,  while  the  BDO  solitons  are  algebraic  Lorentzian  functions.  (8) 
Solution  of  the  TG  equation  is  an  eigenvalue  problem  resulting  in  a  set  of  modes  for  the 
duct  characterized  by  an  eigenvalue  that  is  a  function  of  cq.  As  the  modal  index  is 


29 


increased,  these  modes  have  increasing  numbers  of  vertical  nodes  and  have  decreasing 
phase  speeds.  The  varicose  mode  of  the  Dewan-Picard  model  would  be  called  mode  2  by 
the  oceanographers,  has  a  single  node,  and  values  of  Co  near  60  m/s.  However,  the  NWT 
predicts  a  lower-order  fast  mode,  mode  1,  having  no  nodes  in  the  vertical  and  Co  near  180 
m/s.  This  mode  is  described  as  sinuous  rather  than  varicose,  since  the  parcel 
displacements  above  and  below  the  center  of  the  duct  are  synchronized  rather  than  being 
in  opposition.  Hence,  it  cannot  lead  to  complimentary  airglow  emissions,  regardless  of 
the  duct  location  relative  to  the  airglow  layers.  (9)  Addition  of  a  Burgers  term 
representing  dissipation  to  the  nonlinear  wave  equation  can  produce  turbulent  or  foaming 
bores.  As  the  coefficient  of  the  Burgers  term  is  turned  up,  undulations  disappear  and  the 
frontal  behavior  becomes  monotonic,  while  the  front  still  maintains  its  sharpness. 

Sudden  switching  of  the  strength  of  the  Burgers  term  will  cause  undular  bores  to  evolve 
into  foaming  bores,  or  vice  versa. 


BDO  Field  Snapshots:  T  =  0(10)800 


Figure  13.  Solution  of  the  BDO  equation  for  the  bore  dimensionless  stream  function  as  a  function  of 
dimensionless  horizontal  distance  X  relative  to  a  coordinate  frame  moving  horizontally  with  the  linear 
wave  phase  speed.  The  time  development  is  shown  starting  from  the  initial  value  in  red  as  a  function  of 
the  dimensionless  time  T.  The  time  increment  is  10,  with  the  solutions  for  successive  times  displaced 
upward  by  constant  increments.  The  solutions  are  damped  quasi-sinusoidal  wave  trains  lengthening  with 
time. 

An  example  of  solutions  of  the  BDO  equation  in  this  model  for  small-amplitude 
initial  values  is  shown  in  Figure  13. 
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5.2.  Wave  and  Other  Structures 


5.2.1.  General 

The  AFRL  optical  structure/signatures  team  has  cooperated  with  the  AFOSR/NSF 
Maui  MALT  program,  which  deployed  an  impressive  array  of  sensors  (CCD  imagers, 
spectrometers,  lidar,  radar)  at  the  AMOS  site  for  a  concerted  attempt  to  characterize  the 
sources  of  atmospheric  structure.  Participation  in  prior  campaigns  at  AMOS  (ALOHA 
90,  ALOHA  93)  has  been  very  successful  and  yielded  valuable  insights  into  the 
phenomena  and  mechanisms  that  lead  to  fluctuations  in  optical  emissions.  In  addition, 
work  with  the  MSX  database  has  generated  follow-on  programs  dealing  with  the 
emissions  from  wave-structured  atmospheres. 

Collaborations  with  field  programs,  including  Maui  MALT,  and  outgrowths  of  MSX 
investigations  have  led  to  a  number  of  important  results  relevant  to  the  structure  of  the 
atmosphere  and  the  radiance  it  emits.  Results  on  GWs  are  discussed  in  Section  5.2.2. 

Thunderstorms  are  an  important  source  of  atmospheric  GWs,  and  they  also  generate 
sprites  (upward  lightning  from  cloud  top  to  space);  simultaneous  observations  of  GWs 
and  sprites  are  treated  in  Section  5.2.2.  Section  5.2.3  treats  structuring  by  tides  and 
Section  5.2.4  structuring  by  turbulence. 

5.2.2.  Atmospheric  Gravity  Waves,  Sprites 

The  observation  at  4.3  pm  in  an  MSX  BTH  scene  of  a  circular  GW  pattern  radiating 
from  a  compact  convective  storm  has  been  noted  in  Section3.2.  The  angle  of  propagation 
seen  in  this  pattern  has  been  used  [Dewan  et  al.,  1998]  to  infer  the  wave  frequency  or 
period  and,  with  the  measured  horizontal  wavelength,  to  infer  the  phase  speed  and  the 
transit  time  from  the  source  region  to  the  region  of  observation.  However,  the  MSX  data 
can  be  used  to  infer  further  features  of  the  waves.  In  particular,  the  contrast  of  the 
circular  IR  radiance  pattern  from  unsaturated  GWs  generated  by  an  isolated  thunderstorm 
[Alexander  and  Picard,  unpublished]  depends  on  the  angle  between  the  sensor  LOS  and 
the  GW  phase  front,  maximizing  when  the  LOS  is  aligned  with  the  wavefront.  This  has 
been  called  the  "Venetian-blind"  effect  and  occurs  because  in  the  cases  in  question  the 
CF  is  much  wider  than  the  GW  wavelength  along  the  LOS.  (For  the  MSX  SPIRIT  III 
4.3-mm  B1  band,  the  width  of  the  CF  is  15-20  km,  quite  wide  compared  to  typical 
airglow  layer  widths.)  This  information  can  be  used  to  refine  the  estimate  of  wave 
characteristics  and  explains  certain  features  of  the  observations.  For  example,  it  explains 
why  only  waves  propagating  toward  the  sensor  and  no  waves  propagating  away  from  the 
sensor  are  observed.  This  effect  is  shown  clearly  in  Figure  14. 
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SATELLITE  MOTION 


Modeled  Wave  Radiance 


Figure  14.  Illustration  of  the  way  the  response  of  a  "thick-screen"  IR  emitter  varies 
according  to  the  degree  of  alignment  of  the  LOS  with  the  wavefront.  The  top  panel  shows  a 
thunderstorm-generated  wave  simulation  (courtesy  of  M.  J.  Alexander)  with  several  superposed 
LOS  from  a  satellite-based  sensor.  The  GWs  are  propagating  outward  and  upward  from  the 
center  of  the  storm.  The  bottom  panel  simulates  the  4.3-jnm  VER  from  MSX  as  a  function  of 
nadir  angle  over  the  1°  MSX  FOV,  centered  on  the  angle  where  the  LOS  is  aligned  with  the  phase 
front.  The  simulation  is  run  for  an  unperturbed  atmosphere  and  for  an  atmosphere  perturbed  by 
a  GW.  The  vertical  (horizontal)  wavelength  of  the  wave  is  taken  to  be  10  km  (32  km),  and,  for  a 
sensor  at  the  MSX  orbital  altitude  of  900  km,  perfect  alignment  will  occur  for  a  57°  nadir  angle. 
The  difference  between  the  response  to  the  perturbation  when  the  LOS  is  aligned  (GW 
propagates  toward  sensor)  and  when  it  is  not  aligned  (GW  propagates  away  from  sensor)  is  clear. 

On  the  basis  of  the  MSX  observations  we  received  a  NASA  Office  of  Space  Science 
grant  (PI:  Jeremy  Winick;  period  of  performance:  Apr  2000  -  Mar  2002)  to  develop 
techniques  to  sense  gravity  waves  remotely  using  infrared  emissions  and  airglow  layers 
[* Picard  et  al. ,  2002].  This  work  resulted  in  papers  at  several  conferences  and  a  NASA 
Office  of  Space  Science  report.  One  result  from  this  work  concerns  techniques  for  doing 
rapid  calculations  of  the  atmospheric  MWIR  response  to  GWs.  It  was  demonstrated  that 
an  approach  to  calculating  the  radiative  response  to  GWs  based  on  CFs  for  the 
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unperturbed  atmosphere  works  quite  well  and  cuts  down  enormously  on  computation 
time. 

AFRL  optical  structure/signatures  team  members  Jeremy  Winick  and  Richard  Picard 
have  served  on  the  Waves  Explorer  Science  Team.  Waves  Explorer  is  a  proposed 
$180M  MIDEX  (mid-sized  Explorer)  satellite  (PE  Gary  Swenson,  Univ.  Illinois)  that  will 
carry  out  global  gravity-wave  surveys  and  will  measure  waves  in  the  stratosphere  using 
the  4.3-pm  IR  technique  pioneered  by  MSX.  This  mission  has  been  proposed  to  NASA 
twice  already  (Aug  1998  and  Oct  2001)  without  success.  It  may  be  proposed  again  if  and 
when  the  opportunity  arises,  this  time  with  David  Fritts  (Colorado  Research  Assoc.)  as 
PI. 

CCD  imager  observations  of  715-920  nm  emission  from  Bear  Mtn.,  South  Dakota  in 
support  of  the  NASA  Sprites99  Balloon  Campaign  resulted  in  the  first  simultaneous 
observations  from  the  ground  of  the  sprites  and  GWs  generated  by  a  thunderstorm 
[**Sentman  et  al.,  2003].  The  GWs  were  seen  in  OH  Meinel  high  overtones  while  the 
sprite  emission  resulted  from  the  N2  First-Negative  fundamental  and  overtone  band 
sequences.  In  addition,  this  was  the  first  observation  of  the  complete  circular  concentric 
ring  pattern  of  a  thunderstorm.  We  determined  that  this  was  possible  because  the  GW 
vertical  wavelength  (-20  km)  was  much  greater  than  the  width  of  the  emission  layer  (8 
km).  Thus  there  was  no  Venetian-blind  effect.  In  the  first  observation  of  the  concentric 
ring  pattern  from  thunderstorm-generated  waves  [ Taylor  and Hapgood,  1988]  and  in  the 
MSX  MWIR  observation  from  space  [D ewan  et  al.,  1998],  only  partial  circular  patterns 
were  observed  because  the  vertical  wavelength  was  too  large.  The  total  energy  deposited 
by  the  sprites  in  the  mesosphere  during  the  Sprites99  Campaign  observation  was  -3.3 
MJ,  which  is  well  below  the  level  required  to  have  a  significant  effect  on  the  temperature 
in  the  OH  emission  region  (-1  GJ).  A  video  of  the  GW-sprite  sequence  is  available. 
Several  advances  have  been  made  in  the  description  of  GW  power  spectra.  **Dewan  and 
Grossbard  [2000]  showed  that  published  balloon  data  on  saturated  gravity-wave  power 
spectra,  including  several  recently  published  papers,  were  contaminated  by  spectral 
leakage.  This  resulted  in  enhanced  amplitude  fluctuations  and  a  coupling  between 
spectral  amplitude  and  slope,  with  resulting  implications  for  the  mechanisms  responsible 
for  the  shape  of  saturated  GW  power  spectra. 

The  linear  saturation  theory  proposed  by  Dewan  and  Good  [1986]  to  explain  GW 
spectra  was  tested  in  a  number  of  published  works.  Allen  and  Vincent  [1995]  used 
balloon-borne  temperature  data  to  test  the  model.  They  found  that  the  data  agreed  well 
with  the  work  of  S.  A.  Smith  et  al.  [1987],  who  had  extended  the  Dewan-Good  model,  in 
the  stratosphere,  but  did  not  agree  in  the  troposphere.  Dewan  (unpublished,  2006) 
showed  that  the  lack  of  agreement  occurred  because  S.  A.  Smith  et  al.  [1987]  had  used  a 
separable  power-spectral  model,  unlike  the  unseparable  PSD  in  the  saturated-cascade 
model  of  Dewan  [1997],  When  this  assumption  was  eliminated  and  some  additional 
considerations  were  applied,  the  disagreement  between  theory  and  data  was  eliminated. 

In  addition  to  research  on  second  moments  of  the  saturated  GW  fluctuation  amplitude  (as 
typified  by  GW  PSDs),  important  work  was  also  carried  out  on  properties  of  higher 
moments  of  the  GW  probability  density  function.  A  high-resolution  study  of 
stratospheric  winds  by  smoke  trails  over  scales  where  GWs  are  expected  to  dominate  has 
shown,  for  the  first  time,  that  the  saturated  GW  field  is  intermittent  [* Dewan  and 
Grossbard,  2007].  The  data  is  described  by  Dewan  et  al.  [1984,  1988].  Figure  15 
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compares  the  smoke  trail  data  with  simulated  data  from  a  Gaussian  random  process.  The 
difference  is  clear  with  the  data  showing  large  intermittence,  while  the  Gaussian  process 
is  not  (and  cannot  be)  intermittent.  A  large  "flatness"  of  the  probability  density  function 
is  taken  as  a  quantitative  measure  of  intermittence,  where  flatness  is  defined  as  the  fourth 
central  moment  of  the  distribution  divided  by  its  variance  squared.  It  is  also  shown  that 
the  rms  shears  are  log-normally  distributed.  The  theoretical  interpretation  of  this 
behavior  in  terms  of  Dewan's  saturated-cascade  GW  theory  [Dewan,  1997]  was 
compared  to  that  arising  from  Lindborg's  wave-free  stratified  turbulence  cascade  theory 
[Lindborg,  1999],  While  neither  approach  was  deemed  to  be  completely  satisfactory  yet, 
the  preponderance  of  evidence  appeared  to  favor  Dewan's  saturated-cascade  theory. 


Figure  15.  Illustration  of  intermittence  in  high-resolution  (10-m)  wind  profile  data  (left  panel).  Data 
shown  are  vertical  shears  of  total  horizontal  wind  (s'1),  plotted  from  15  to  45  km  altitude.  The  right  panel 
shows  the  shear  for  a  simulated  Gaussian  random  process  wind  profile  for  comparison.  Both  sets  of  shears 
have  similar  power  spectra,  but  the  tendency  of  the  real  data  to  stay  near  zero  most  of  the  time  with 
sporadic  very  large  excursions  contrasts  with  the  small  but  continuous  excursions  of  the  simulated  data. 


5.2.3.  Tides 

The  ARC  non-LTE  radiative-transfer  model  has  been  modified  to  incorporate  tidal 
modulation  of  the  atmospheric  temperature  and  density  in  the  calculation  of  the  15-pm 
CO2  V2  emission.  This  is  necessary  in  order  to  model  reliably  the  LT  dependence  of 
observed  radiance  structure.  See  the  discussion  in  Sections  4.3  and  4.5  describing  the 
occurrence  of  thermospheric  radiance  layers  or  “knees”  in  the  SABER  data  and  the 
comparison  of  the  data  to  this  model. 


5.2.4.  Turbulence 

The  role  of  turbulence  in  structuring  light  waves  propagating  in  the  free  atmosphere 
has  been  addressed  by  Dewan  et  al.  [1993],  The  so-called  Dewan  model  predicts  the 
altitude  profile  of  the  optical  refractive-index  structure  constant  Cn2,  which  depends  on 
high-spatial-resolution  temperature  measurements,  on  the  basis  of  radiosonde  data  only, 
which  have  inherently  low  spatial  resolution.  This  model  has  important  application  to 
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determining  ground-based  astronomical  seeing  [*Dewan,  2003]  and  has  resulted  in  an 
invention  disclosure  by  Dewan  et  al.  For  this  achievement,  team  member  Edmond 
Dewan  received  the  2002  Harold  Brown  Award  from  the  Secretary  of  the  Air  Force  in 
May  2003. 

In  a  recent  publication  [ **Dewan  and  Grossbard,  2007],  it  has  been  shown  that  a 
result  of  Tatar  ski  [1971]  relating  Cn  to  the  outer  scale  of  turbulence  Lo  in  the 
atmospheric  boundary  layer  (ABL)  severely  underpredicts  Cn  in  the  free  atmosphere  and 
must  be  modified.  The  reason  for  the  underprediction  involves  the  smaller-than-expected 
thickness  of  the  turbulent  layers  and  the  minor  fraction  of  the  free  atmosphere  that  is 
turbulent. 

Turbulence  fields  are  clearly  intermittent.  Recall  also  the  results  of  Section  5.2.2, 
where  GWs  were  also  shown  to  be  intermittent.  Suppose  it  should  happen  that 
Lindborg's  theory  [Lindborg,  1999]  cited  there  that  interacting  saturated  GWs  form  a 
stratified  turbulence  cascade  is  correct.  Then  the  range  of  spatial  scales  that  can  be 
considered  turbulent  will  be  extended  considerably  at  the  large-scale  end,  and  the  results 
on  GW  intermittence  cited  above  could  be  considered  results  on  the  intermittence  of 
(extended-range)  turbulence. 


6.  RADIATIVE  TRANSFER  AND  REMOTE  SENSING 
6.1.  Radiative  Transfer  Methods 

A  detailed  intercomparison  of  three  different  world-class  4. 3 -pm  non-LTE  radiative 
models  and  an  operational  non-LTE  model  has  been  carried  out.  These  codes  are  the 
AFRL  ARC  line-by-line  (LBL)  radiative-transfer  model  [  Wintersteiner  et  al.,  1992; 

Nebel  et  al.,  1994],  the  IAA  (Institute  de  Astro fisica  de  Andalucia)  Curtis-matrix  narrow- 
band  model  [Lopez-Puertas  et  al.,  1986;  Lopez-Puertas  and  Taylor,  2001],  the  hybrid 
SABER  research  model  based  on  LINEPAK  [ Gordley  et  al.,  1994]  and  developed  by 
GATS  (G  &  A  Technical  Software)  and  NASA  Langley  Research  Center  (LaRC),  and 
the  GATS/LaRC  SABER  operational  non-LTE  radiative  code,  based  on  BANDPAK 
[. Marshall  et  al.,  1994],  The  agreement  between  the  results  obtained  with  these  codes 
was  remarkable,  and  any  sources  of  discrepancy  were  isolated  and  identified.  Since  these 
codes  have  different  heritage,  were  developed  independently,  and  use  independent  and 
entirely  different  methods/algorithms,  the  agreement  between  them  serves  as  evidence 
that  the  processes  involved  in  this  complex  emission  are  well  understood.  This  has 
resulted  in  confidence  in  the  AF  radiative  codes  and  their  ability  to  predict  MWIR  (mid¬ 
wave  infrared)  backgrounds.  In  earlier  work  [Lopez-Puertas  et  al.,  1994]  the  ARC  and 
IAA  non-LTE  15-pm  models  were  compared,  with  similar  good  results. 

LBL  radiative-transfer  methods  are  the  most  accurate  methods  for  radiative  transfer 
calculations.  They  are  devoid  of  simplifying  assumptions  about  the  nature  of  the 
emission  and  free  of  the  many  pitfalls  affecting  band-model  approaches.  However,  LBL 
calculations  are  notoriously  time  consuming.  Methods  for  increasing  the  speed  of  non- 
LTE  radiative  models  and  generating  rapid  LBL  radiative-transfer  algorithms  have  been 
developed  and  implemented  in  the  ARC  model.  These  methods  make  possible  otherwise 
unfeasible  LBL  calculations,  enabling  very  rapid  and  accurate  computation  of  vibrational 
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populations  for  different  conditions.  These  acceleration  methods  depend  on  the  use  of 
stored  functions  to  facilitate  calculations  and  result  in  improved  convergence  and 
accuracy  of  the  lambda-iteration  technique  used.  We  also  found  it  important  to  correct 
properly  for  the  assumption  of  homogeneous  layers  when  considering  optically  thick 
bands.  Using  CO2  emissions  as  an  example,  we  have  computed  rapidly  obtained  radiance 
and  cooling  results  that  are  useful  to  characterize  the  emissions  from  evolving  and 
variable  scenes.  The  algorithms  were  extended  to  time- varying  situations,  exemplified 
by  auroras  and  have  also  been  adapted  to  atmospheres  perturbed  by  GWs,  as  discussed  in 
Section  5.2.2. 

Many  improvements  in  the  physics  and  photochemistry  have  been  made  in  the  ARC 
non-LTE  LBL  model  for  atmospheric  CO2  4.3-pm  and  1 5 -pm  emission,  including 
addition  of  several  V-V  transfer  processes.  Note  also  the  discussion  in  Section  4.2  of  the 
addition  of  V-V  coupling  terms  redistributing  energy  among  CO2  isotopes.  As  discussed 
there,  addition  of  this  process  solved  the  problem  of  the  discrepancies  between  SABER 
temperature  retrievals  in  the  summer  polar  mesopause  region  and  data,  and  this  change 
was  kept  as  a  permanent  modification  to  ARC. 

Finally,  team  members  participated  as  consultants  and  evaluators  in  the  formulation 
of  the  first  book  devoted  to  non-LTE  radiative  transfer  in  atmospheres.  The  book  is 
entitled  Non-LTE  Radiative  Transfer  in  the  Atmosphere,  and  its  authors  are  Manuel 
Lopez-Puertas  of  IAA,  Granada,  Spain  and  Frederick  W.  Taylor  of  Oxford  University. 
The  book  discusses  AFRL  models,  in  addition  to  treating  AF  and  DoD  space-based 
measurement  programs  in  detail,  starting  with  ICECAP  [ Stair  et  al.,  1975]  and  SPIRE 
[Stair  et  al.,  1985]  and  continuing  through  CIRRIS-1A  and  MSX.  Participation  by  team 
members  ensured  that  AF/DoD  programs  were  properly  represented  and  helped  to  assure 
the  quality  of  the  work. 


6.2.  Radiative  Transfer  in  Structured  Atmospheres 

Our  research  in  radiative  transfer  in  structured  atmospheres  is  discussed  in  Sections 
3.2,  3.3,  4.5,  5.2.3,  and  6.1.  A  review  of  some  of  this  work  has  been  published  [* Picard 
et  al.,  2002], 


6.3.  Temperature  and  Species  Profile  Retrieval  from  Remote 
Sensing  Data 

In  the  upper  stratosphere  and  lower  mesosphere,  in  daytime  the  CO2  laser  bands  near 
10  pm  (transition  0001 1  -  10001  in  HITRAN  [Rothman  et  al.,  2005]  or  Benedict 
notation)  are  populated  by  both  LTE  processes  (thermal  collisions)  and  non-LTE 
processes  [solar  fluorescent  pumping,  pumping  by  upwelling  earthshine,  energy  transfer 
from  N2(v)  excited  by  0(*D)].  The  high-resolution  (~1  cm'1)  IR  limb  emission  measured 
by  the  DoD  CIRRIS-1 A  spectrometer  flying  on  the  Space  Shuttle  STS-39  clearly 
illustrates  these  effects.  The  presence  of  non-LTE  emissions,  if  left  unaccounted  for, 
would  seriously  compromise  the  use  of  the  laser-band  emission  to  retrieve  kinetic 
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temperature  profiles.  Since  ozone  radiates  in  this  same  spectral  region,  non-LTE 
emissions  could  also  compromise  the  retrieval  of  ozone  vmr  from  the  data,  particularly  in 
the  daytime.  A  non-LTE  compensation  technique  was  applied  to  the  data  whereby  the 
non-LTE  component  of  the  radiance  was  modeled  and  subtracted  from  the  total  measured 
radiance  [**Miller  et  al.,  2000],  The  kinetic-temperature  profile  and  the  O3  vmr  were 
then  derived  using  the  modified  data.  The  correction  for  non-LTE  used  the  ARC  model 
[Nebel  et  al.,  1994]  to  calculate  the  CO2  Tv;b  and  FASCODE  [Smith  et  al.,  1978; 
Anderson  et  al.,  1994]  to  compute  the  LOS  radiative  transfer.  The  LTE  retrievals  were 
then  carried  out  on  the  residual  radiance  by  Rodgers'  optimal-estimation  method 
[Rodgers,  1992;  2000]  using  FASCODE  as  the  forward  LTE-radiance  model  [Miller  et 
al.,  1999].  The  bands  used  in  the  retrievals  were  940-960  cm'1  (CO2  laser-band  P  branch) 
and  789-794  cm'1  (first  laser  hot  band  Q  branch)  for  the  temperature  and  1073-1097  cm'1 
(O3  V3)  for  the  ozone  vmr. 

The  SABER  retrieval  problem  was  alluded  to  earlier  in  Section  4.1.  Every  SABER 
channel  possesses  non-LTE  radiators,  either  as  the  target  species  or  as  interfering  species 
or  both.  A  more  fundamental  approach  to  retrieval  was  necessary  in  this  case,  and  the 
SABER  retrievals,  operational  and  otherwise,  use  a  full  non-LTE  model  for  vibrational 
temperatures  as  an  integral  part  of  the  forward  radiative  model.  The  full  non-LTE 
forward  model  is  used  in  a  Rodgers-type  optimal-estimation  retrieval  to  do  the  inversion, 
similarly  to  the  approach  in  the  previous  paragraph. 

6.4.  Foliage-Canopy  Radiative  Transfer 

An  atmospheric  radiative-transfer  code  (MODTRAN  4)  has  been  interfaced  to  a 
foliage-canopy  model  (LCM2)  for  the  first  time,  enabling  more  accurate  predictions  of 
terrestrial  sub-canopy  structures  seen  from  space  [* Picard  et  al.,  2001;  Ganapol  et  al., 
*2002,  *2003].  The  optical/infrared  radiation  emanating  from  the  top  of  a  foliage  canopy 
contains  spectral  information,  and  perhaps  also  polarization  and  temporal  information, 
about  targets  and  human  activity  within  and  beneath  the  canopy,  even  when  the  foliage  is 
totally  obscuring.  We  have  undertaken  a  modeling  study  to  determine  the  extent  to 
which  this  information  may  be  recovered  from  hyperspectral  imagery  acquired  by  space- 
based  sensors.  Such  passive  methods  hold  the  promise  of  considerable  size,  weight,  cost, 
and  power  advantages  over  foliage -penetration  radar.  Using  methods  from  radiative- 
transfer  theory,  we  have  found  that  the  spectrally  resolved  upwelling  radiance  can 
indicate  the  presence  and  motion  of  a  simulated  target  beneath  a  canopy  under  totally 
obscuring  foliage,  that  is  where  the  optical  depth  of  the  canopy  is  large  or  the  leaf  area 
index  (LAI)  ~5. 

Using  the  LCM2  canopy  code  developed  by  collaborator  Barry  Ganapol  of  Univ. 
Arizona  [Ganapol  et  al.,  1999],  we  have  emphasized  rapid  exact  analytic  solutions  of  the 
equation  of  transfer  based  on  the  Fn  (Facilitation-N)  method  [Siewert  et  al.,  1980].  More 
intensive  discrete-ordinates  or  Sn  methods  (for  example,  Evans  and  Marshak  [2005]) 
were  also  used  along  with  computation  acceleration  techniques.  LCM2  describes 
radiative  transfer  at  two  nested  levels— within  an  individual  leaf  (LEAFMOD)  as  well  as 
within  a  cloud  of  oriented  leaves  forming  the  canopy  (CANMOD).  LEAFMOD  makes 
use  of  the  specified  absorption  and  scattering  properties  of  leaves,  while  CANMOD  uses 
the  LAI  and  leaf-angle  distribution  (LAD)  of  the  canopy  as  well  as  the  output  of 
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LEAFMOD.  The  lower  boundary  beneath  the  canopy  is  an  absorptive/reflective  surface 
simulating  soil  or  a  man-made  target  material. 

The  AF  standard  atmospheric  radiative  band  model  MODTRAN  was  used  to  describe 
the  radiative  transfer  from  the  top  of  the  canopy  (TOC)  to  the  satellite  [ Berk  et  al.,  1999, 
2004],  We  have  focused  on  the  solar  reflective  regime  (0.4  -  2.5  pm),  including  the 
visible  (VIS),  near-infrared  (NIR),  and  short-wave  IR  (SWIR),  and  have  found 
considerable  differences  between  the  behavior  in  the  VIS  and  the  NIR  range.  Photons 
strike  the  TOC,  are  multiply  scattered  and  absorbed  in  canopy,  reflect  from  surfaces 
beneath  canopy,  traverse  the  canopy  again  with  multiple  scattering  and  absorption,  and 
reemerge  from  the  TOC.  In  so  doing,  they  carry  the  spectral  signature  of  elements  in  the 
canopy  and  targets  beneath  it.  In  the  VIS  region,  the  canopy  is  strongly  absorptive  and 
target  signatures  are  weak.  However,  beyond  the  chlorophyll  absorption  edge,  in  the  NIR 
and  the  SWIR,  the  absorption  decreases  dramatically,  and  a  relatively  large  fraction  of 
the  photons  incident  on  the  canopy  reemerge  from  the  TOC  and  propagate  to  the 
spacecraft  sensor  imprinted  with  the  signature  of  the  target  materials.  This  is  true 
provided  one  avoids  spectral  regions  where  H20  absorption  bands  of  the  canopy  and 
atmosphere  lie.  The  same  is  true  in  certain  regions  of  the  SWIR,  once  again  avoiding 
water  bands,  although  the  incident  solar  radiance  is  declining  fast  here. 

Under  the  SFFP  (Summer  Faculty  Fellowship  Program),  LCM2  has  now  been 
extended  to  include  a  treatment  of  polarization  [ *Ganapol  et  al.,  2005].  Specular 
reflection  has  been  introduced  into  the  leaf  scattering  phase  function  in  LCM2. 

Laboratory  and  field  studies  have  demonstrated  that  the  specular  component  of  scatter 
from  the  leaf  surface  is  the  primary  source  of  linear  polarization  in  a  canopy  signal 
[] Vanderbilt  et  al.,  1985a,  1985b].  So,  first  linear  polarization  was  introduced,  followed 
by  progress  toward  a  full  vector  theory  of  polarization. 


6.5.  Remote  Sensing  Meeting  Organization 

Team  member  Picard  is  principal  Chair  of  the  Conference  on  Remote  Sensing  of 
Clouds  and  the  Atmosphere  to  be  held  as  part  of  the  SPIE  Europe  Remote  Sensing  and 
Security  &  Defence  symposia  in  Cardiff,  UK  in  Sep  2008.  He  has  served  as  a  Co-Chair 
and  Program  Committee  member  for  the  conference  since  2001,  introducing  for  the  first 
time  a  focus  on  remote  sensing  of  the  middle  and  upper  atmosphere  to  the  conference. 
Other  topics  represented  at  the  conference  include  radiative  transfer  and  various  other 
aspects  of  ground-based  and  space-based  remote  sensing.  See,  for  example,  Comeron  et 
al.  [2007]  for  the  proceedings  of  the  2007  conference  in  Florence,  Italy.  Previous 
conferences  were  held  in  Brugge,  Belgium  (2006);  Stockholm  (2005);  Gran  Canaria, 
Spain  (2004);  Barcelona  (2003);  Agia  Pelagia,  Crete  (2002);  and  Toulouse,  France 
(2001).  See  the  proceedings  at  Schafer  et  al.  [*2003,  *2004a,  *2004b,  *2005], 
Participation  in  organizing  the  conference  enables  AF  access  to  the  international  remote¬ 
sensing  and  structured-atmosphere  research  communities.  The  symposia  address  issues 
in  the  DoD  and  foreign  defense  communities,  as  well  as  NASA/ESA  issues. 
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7.  SPACE  SURVEILLANCE 


In  the  later  years  task  activities  have  been  expanded  and  reoriented  in  part  to  support 
remote  sensing  and  signature  studies  of  interest  for  space  situational  awareness  (SSA) 
and  space-object  identification  (SOI)  applications.  The  purpose  of  this  activity  is  to 
investigate  the  spatial,  spectral,  and  temporal  signatures  characterizing  RSOs. 

We  have  considered  space  objects  illuminated  by  non-conventional  sources,  including 
upwelling  earth  and  earth-atmosphere  radiation  (earthshine  and  airglow).  In  addition  to 
influencing  the  environment  of  RSOs  and  providing  a  source  of  radiance  from  which  they 
must  be  distinguished,  airglow  and  earthshine  can  illuminate  RSOs  in  eclipse  or  during 
the  dimmer  lunar  phases  and  at  new  moon.  OH  Meinel-band  airglow  is  attractive  as  an 
illuminant  because  its  radiance  when  integrated  across  all  spectral  bands  is  many  mega- 
Rayleighs.  See,  for  example,  Gruninger  et  al.  [2007].  Since  its  strongest  bands  are  in  the 
vibrational  first-overtone  region  near  1.5  pm,  a  target  region  for  optical  fiber 
communications  because  of  the  fiber  extinction  minimum  in  that  region,  one  should  be 
able  to  take  advantage  of  advances  in  NIR  detector  technology  to  make  use  of  OH  Meinel 
airglow  more  attractive.  This  approach  is  particularly  appealing  for  objects  in  low  earth 
orbit  (LEO). 

The  OH  airglow  source  has  been  implemented  in  the  Optical  Signatures  Code  (OSC). 
Although  OSC  was  developed  over  an  extended  period  by  Teledyne  Brown,  the  Army 
Missile  Command,  and  predecessor  organizations  for  missile-defense  applications 
r www.teledvnesolutions.com/docs/specsheets/041 4  Optical  Siqnatures.pdfl.  it  contains  a  lot 
of  the  functionality  required  for  simulating  the  radiative  signatures  of  space  objects. 

Because  of  the  long  optical  path  through  a  layer  in  the  limb,  the  limb  emission  is  most 
important  in  determining  the  RSO  illumination.  However,  the  limb  enhancement  peak  is 
very  narrow.  On  the  other  hand,  the  sub-limb  (BTH)  airglow  comes  from  shorter  paths 
but  is  associated  with  larger  solid  angles.  An  accurate  model  must  take  this  subtlety  into 
account.  The  full  description  of  the  problem  is  subject  to  the  intrinsic  diurnal,  seasonal, 
latitudinal,  and  solar-cycle  variability  of  the  OH  airglow,  which  is  catalogued  by  SABER. 
The  SABER  database  contains  an  extensive  multi-year  catalog  of  volume  emission  rates 
in  two  bright  SWIR  OH  Meinel  bands  of  interest  for  this  application.  As  a  result  of  the 
variability  of  the  OH  airglow,  there  can  be  significant  variation  of  the  upwelling  radiance 
impinging  on  the  RSO.  One  would  expect  the  impinging  radiance  to  vary  as  the  RSO 
and  the  sub-RSO  point  move.  This  variation  is  due  to  a  mixture  of  photochemical  effects 
plus  dynamical  (tidal  and  other)  effects.  So  a  parametrized  climatology  of  this  variation 
from  the  SABER  database  may  be  of  value. 

When  a  ground-based  sensor  rather  than  a  space  sensor  is  viewing  an  eclipsed  RSO, 
the  use  of  airglow  as  an  illuminant  is  much  less  attractive.  This  is  because  one  must  look 
through  the  airglow  layer  in  the  foreground,  which  will  reduce  the  contrast  of  the  RSO- 
reflected  airglow  against  the  background/foreground  airglow. 

Other  airglow  bands  and  earthshine  bands  from  CO2,  H2O,  NO,  and  other  species 
radiating  at  wavelengths  less  than  2.5-3  pm  can  also  be  considered  as  illuminants.  At 
longer  wavelengths  in  the  thermal  regime,  however,  the  atmosphere  and  the  RSO  become 
self-emitting,  and  an  illuminant  is  so  not  useful. 

We  have  also  started  to  consider  methods  of  RSO  probing  by  using  non-conventional 
methods.  The  latter  include  exploiting  the  information  found  in  multipoint  and  multi- 
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time  correlation  functions.  This  leads  to  the  consideration  of  interferometric  methods 
(for  example,  Monnier  [2003]),  including  both  amplitude  and  intensity  interferometry,  for 
resolved  and  unresolved  objects.  Temporal  interferometers,  for  example  the  prototype 
Michelson  two-beam  interferometer,  have  the  potential  to  detect  quite  subtle  temporal 
changes  in  targets,  including  those  due  to  rotation,  vibration,  and  relative  motions. 

Spatial  interferometers,  such  as  the  prototype  Michelson  stellar  interferometer,  can 
increase  spatial  resolution  by  effectively  forming  synthetic  apertures.  In  addition,  we 
have  turned  our  attention  more  general  techniques  to  enhance  spatial  resolution.  Called 
generally  superresolution  techniques  (for  example,  Zalevsky  and  Mendlovic  [2004]), 
these  methods  seek  to  beat  the  so-called  "diffraction  limit"  or  the  limitations  imposed  by 
the  usual  optical- system  point-spread  function. 


8.  COLLABORATIONS 
8.1.  General 

Active  collaborations  have  occurred  with  NASA  (TIMED/SABER  satellite,  Waves 
Explorer  satellite,  MSX  data  exploitation  and  gravity-wave  remote  sensing)  and  with 
MDA  (MSX  modeling  and  analysis,  radiative  code  enhancements  for  SBIRS). 
Collaborations  have  also  occurred  with  NRO  and  the  AF  Warfighter  hyperspectral 
imaging  satellite  program  (foliage  canopy  and  atmospheric  radiative  transfer),  with  the 
AFOSR/NSF  Maui  MALT  program  (gravity  waves  and  localized  structures),  and  with 
ESA  (non-LTE  radiative  codes,  ENVISAT  1  /  MIPAS).  These  collaborations  are 
generally  treated  in  the  preceding  sections.  This  section  ends  with  a  discussion  of  a  few 
collaborations  not  addressed  above. 


8.2.  Support  to  Columbia  Shuttle  Accident  Investigation 

Team  members  Richard  Picard  and  Edmond  Dewan  contributed  to  the  investigation 
following  the  Columbia  Space  Shuttle  accident  in  Feb  2003  on  behalf  of  AFRL.  Picard 
served  on  a  distinguished  panel  of  university  and  government  experts  investigating 
space/atmospheric  environment  effects  on  the  Columbia  accident  and  coordinated 
drafting  of  the  panel’s  report  [* Picard  et  al.,  2003].  The  panel  was  convened  at  the 
request  of  the  NASA  Johnson  Space  Center’s  Space  Shuttle  Vehicle  Engineering  Office, 
and  its  report  concludes  that  an  image  that  claimed  to  show  a  lightning  bolt  striking  the 
reentering  Columbia  orbiter  was  the  result  of  an  artifact  and  that  space/atmospheric 
environment  effects  by  themselves  were  not  responsible  for  the  accident.  Dewan 
provided  a  tutorial  on  weather  in  the  mesospheric  region  of  the  atmosphere  (50-85  km 
altitude)  to  the  Columbia  Accident  Investigation  Board  at  their  request.  This  region  can 
be  subject  to  large  temperature  gradients  and  large  wind  shears,  as  well  as  ambient 
atmospheric  electrical  effects,  and  is  critical  to  the  reentering  Shuttle  orbiter,  in  this  case 
representing  the  region  where  the  orbiter  broke  up  on  reentry. 
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8.3.  High  Altitude  Airship  (HAA)  and  Near-Space  Study 


The  potential  hazards  for  high-altitude  airships,  such  as  envisioned  in  MDA’s  HAA 
program,  and  their  systems  due  to  upward  lightning  from  cloud  tops  were  assessed.  The 
little-studied  streamer  discharges  such  as  blue  jets  [  Wescott  et  al.,  1995,  1998],  blue 
starters  [Wescott  et  al.,  1996],  and  the  recently  discovered  gigantic  jets  [Pasko,  2003;  Su 
et  al.,  2003]  are  of  particular  interest  as  potential  hazards.  In  addition,  the  potential  for 
using  such  airships  as  platforms  for  investigating  the  emission  from  high-altitude 
discharges  and  other  background  sources  was  considered.  Support  was  also  provided  to 
the  AFSPC/AFRL  Near-Space  Study. 


8.4.  C/NOFS  and  Triggering  of  Spread-F 

Exploratory  studies  were  carried  out  relating  to  GW  triggering  of  equatorial  spread-F 
(ESF)  by  seeding  the  Rayleigh-Taylor  instability  [Hysell  et  al.,  1990;  Huang  et  al.,  1993; 
Huang  and  Kelley,  1996a,  1996b,  1996c],  The  occurrence  statistics  of  spread-F  suggests 
that  it  could  be  caused  by  GWs  generated  by  convective  storms  [McClure  et  al.,  1998]. 
Since  200-km-scale  structure  is  seen  in  radar  backscatter  [Hysell  et  al.,  1990;  Huang  et 
al.,  1993]  and  in  the  630-nm  airglow  emitted  by  atomic  oxygen  from  the  perturbed  F 
region  thermosphere,  this  suggests  that  the  horizontal  wavelength  Ah  of  any  triggering 
GWs  should  be  -200  km.  Presuming  that  the  triggering  GW  is  generated  in  the 
troposphere  and  propagates  linearly  to  the  equatorial  F  region,  then  Ah  for  the  generated 
wave  must  also  be  -200  km.  For  the  wave  to  propagate  to  the  F  region  without  severe 
attenuation,  the  vertical  wavelength  Az  must  be  large  (-60  km).  It  follows  from  the  GW 
dispersion  relation  that  the  period  x  of  the  wave  must  be  short  (-15  min).  Such  waves  do 
not  appear  to  be  generated  by  typical  squall  lines  [Alexander  et  al.,  1995],  but  they  can  be 
generated  by  intense  large  long-lived  thunderstorm  complexes  with  cold  high  cloud  tops 
(mesoscale  convective  complexes,  or  MCCs)  that  sometimes  develop  from  squall  lines 
[  V a  das  and  Fritts,  2004], 

As  an  alternative,  if  the  Prakash  mechanism  [Prakash,  1999]  is  active  then  perhaps 
the  triggering  GW  only  needs  to  propagate  to  the  extratropical  E  region,  where  it  can 
induce  a  plasma  excitation.  The  plasma  wave  then  maps  the  excitation  along  the  field 
line  to  the  equatorial  F  region,  triggering  ESF.  Waves  having  shorter  Az  can  easily 
propagate  to  the  E  region.  Thus,  with  this  mechanism,  waves  generated  by  typical  squall 
lines  as  well  as  by  MCCs  can  both  propagate  to  the  E  region  and  play  a  role  in  triggering 
spread-F.  Notice,  however,  that  even  though  GWs  generally  propagate  obliquely 
upward,  it  is  still  true  that  one  would  expect  a  larger  horizontal  separation  between  the 
GW  source  region  and  the  location  of  the  spread-F  in  the  case  of  the  Prakash  mechanism. 
Hence,  the  environment  in  which  the  GWs  are  generated  and  propagate  can  be  quite 
different  for  the  two  mechanisms. 

A  procedure  was  also  suggested  for  the  TIMED  spacecraft  to  investigate  the  influence 
of  atmospheric  GWs  on  ESF.  The  procedure  makes  use  of  data  from  two  instruments  on 
TIMED,  specifically,  SABER  and  GUVI.  SABER  data  is  used  to  look  at  GWs  from  15  - 
100  km  altitude,  as  discussed  in  Section  4.9.  The  GUVI  uv  instrument  can  see  spread-F 
on  a  global  basis  using  the  atomic-oxygen  135.6-nm  emission.  In  the  F  region  this  signal 
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is  due  to  radiative  recombination  of  0+  ions  and  thus  is  proportional  to  the  electron 
density  squared. 

These  investigations  make  use  of  team  members'  expertise  in  atmospheric  GWs  and 
support  the  goals  of  the  C/NOFS  program.  It  is  expected  that  they  will  continue  in  the 
future. 


8.5.  Satellite  Drag 

SABER  findings  on  thermosphere/mesosphere  variability  and  structure  were 
reviewed  for  the  benefit  of  the  satellite-drag  research  community.  SABER's  ability  to 
retrieve  remotely  temperature  and  pressure  profiles  as  a  function  of  altitude  on  a  global 
basis  means  that  density  profiles  as  a  function  of  altitude  are  also  known  globally. 

Density  is  a  critical  parameter  determining  drag.  These  retrievals  are  routine  through  the 
mesopause  and  up  to  100  km.  Further  extensions  into  the  lower  thermosphere  up  to  the 
noise  limit  are  also  carried  out,  but  are  dependent  on  the  proper  retrieval  or  modeling  of 
atomic  oxygen. 


9.  CONCLUSION 

This  report  describes  a  varied,  rich,  and  productive  research  program  carried  out  over 
an  extended  period  in  support  of  AF  Space  Situational  Awareness  (SSA)  and  Space- 
Object  Identification  (SOI)  functions.  Results  have  included  significant  advances  based 
on  the  exploitation  of  space-based  remotely  sensed  data,  especially  data  from  MSX  and 
from  TIMED/SABER.  We  have  also  extended  the  knowledge  and  utility  of  non-LTE 
radiative  transfer,  atmospheric  GWs,  inversion  layers,  bores,  and  turbulence.  These 
advances  have  been  recognized  by  the  AF  in  the  form  of  significant  AF  achievement 
awards.  These  advances  indicate  bright  prospects  for  future  impacts  of  related  work  on 
SSA  and  space-object  detection/characterization. 
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AirBome  Laser 

AF 
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Air  Force  Base 

AFOSR 

Air  Force  Office  of  Scientific  Research 
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Air  Force  Research  Laboratory 

ALOHA 

Airborne  Lidar  and  Observations  of  the  Hawaiian  Airglow 
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Air  Force  Maui  Optical  and  Supercomputing  Site 

ARC 

Atmospheric  Radiance  Code 
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Advanced  SPace  ENvironment  (version  of  TIME-GCM  code) 
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Above-the-horizon 

BANDPAK 

(Subroutine  library  for  broadband  atmospheric  radiative  transfer) 
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Benjamin-Davis-Ono  (nonlinear  wave  equation) 

BTH 

Below-the-horizon 

BV 

Brunt-Vaisala  (BV  frequency  =  frequency  of  buoyancy  oscillations) 
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Cornell  All-sky  Roving  Imager 
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Community  Aerosol  and  Radiation  Model  for  Atmospheres 
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Cornell  All-Sky  Imager 

CCD 

Charge-coupled  device 

CF 

Contribution  function 

CIRRIS 

Cryogenic  InfraRed  Radiance  Instrumentation  for  Shuttle 

C/NOFS 

Communication  /  Navigation  Outage  Forecasting  System 

C02 

Carbon  dioxide 

CRISTA 

Cryogenic  Infrared  Spectrometers  and  Telescopes  for  the  Atmosphere 

CSU 

Colorado  State  University 

DNS 

Direct  numerical  simulation 

DoD 

Department  of  Defense 

DT 

Diurnal  tide 

EOS 

Earth  Observing  System 

ESA 

European  Space  Agency 

ESF 

Equatorial  Spread-F 

ENVISAT 

ENVIronmental  SATellite 

FASCODE 

Fast  Atmospheric  Signature  CODE 

FAST 

Fast  Auroral  SnapshoT  Explorer 

FLIP 

Field-Line  Interhemispheric  Plasma  (model) 

FOV 

Field-of-view 

GATS 

G&A  Technical  Software 

GSWM 

Global  Scale  Wave  Model 
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GUVI 

GW 

HAA 

HALOE 

HITRAN 

IAA 

IBC 

ICECAP 

IR 

IRI 

JPL 

KdV 

Kp 


LaRC 

LBH 

LBL 

LCM 

LINEPAK 

LOS 

LT 

LTE 

LWIR 

L0 

MaCWAVE 

MALT 

MCC 

MIDEX 

MAPSTAR 

MDA 

MIPAS 

MLS 

MLT 

MODTRAN 

MSIS 

MSSS 

MSX 

MWIR 


Global  Ultraviolet  Imager 
Gravity  (or  buoyancy)  wave 

High- Altitude  Airship 
HALogen  Occultation  Experiment 

High-resolution  TRANsmission  molecular  absorption  database 
Instituto  de  Astrofisica  de  Andalucia 

International  Brightness  Coefficient  -  index  for  visible  aurora 
Infrared  Chemistry  Experiments  -  Coordinated  Auroral  Program 
Infrared 

International  Reference  Ionosphere 

Jet  Propulsion  Laboratory 

Korteweg-deVries  (nonlinear  wave  equation) 

Planetary  K-index  (3 -hour,  measure  of  horizontal  magnetometer 
fluctuations) 

Langley  Research  Center 
Lyman-Birge-Hopfield  (band  system) 

Line-by-line 
Leaf  Canopy  Model 

(Subroutine  library  for  high-resolution  atmospheric  radiative  transfer) 

Line-of-sight 

Local  time 

Local  thermodynamic  equilibrium 
Long-wave  infrared 
Outer  scale 

Mountain  and  Convective  Wave  Ascending  VErtically 
Mesosphere  and  Lower  Thermosphere  (program,  as  in  Maui-MALT) 
Mesoscale  Convective  Complex 
MID-size  EXplorer 

Middle  Atmosphere  Periodic  Structure  Associated  Radiance 
Missile  Defense  Agency 

Michelson  Interferometer  for  Passive  Atmospheric  Sensing 
Microwave  Limb  Sounder 
Mesosphere  and  Lower  Thermosphere 
MODerate-resolution  TRANsmittance  code 
Mass  Spectrometer,  Incoherent  Scatter  radar  model 
Maui  Space  Surveillance  Station 
Midcourse  Space  Experiment 
Mid-wave  infrared 
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NASA 

National  Aeronautics  and  Space  Administration 

NE 

Physics  and  Electronics  Directorate,  AFOSR 

NIR 

Near  IR 

NO 

Nitric  oxide 

NRLMSISE 

Naval  Research  Laboratory  Mass  Spectrometer,  Incoherent  Scatter 
radar,  Extended  model 

NSF 

National  Science  Foundation 

NWT 

Nonlinear  wave  theory 

OH 

Hydroxyl  radical 

01 

Atomic  oxygen  (neutral) 

OSC 

Optical  Signatures  Code 

OSIRIS 

Optical  Spectrograph  and  InfraRed  Imaging  System 

PI 

Principal  investigator 

PSD 

Power  spectral  density 

Ri 

Richardson  number 

RSO 

Resident  space  object 

RVB 

AFRL,  Space  Vehicles  Directorate,  Battlespace  Environment  Division 

SABER 

Sounding  of  the  Atmosphere  using  Broadband  Emission  Radiometry 

SAG 

SHARC  Atmosphere  Generator 

SAMM 

SHARC  And  MODTRAN  Merged 

SDT 

Semidiurnal  tide 

SFFP 

Summer  Faculty  Fellowship  Program 

SHARC 

Strategic  High- Altitude  Radiance  Code 

SOAR 

Structured  Optical  Atmospheric  Radiance 

SOI 

Space  object  Identification 

SPIM 

SPectral  IMager 

SPIRE 

SPectral  Infrared  Rocket  Experiment 

SPIRIT 

SPectral  (or  SPace)  InfraRed  Imaging  Telescope, 

SSA 

Space  Situational  Awareness 

SSGM 

Strategic  Scene  Generator  Model 

SWIR 

Short-Wave  IR 

TG 

Taylor-Goldstein  (equation) 

TI 

Temperature  inversion 

TIDI 

TIMED  Imaging  Doppler  Interferometer 

TIMED 

Thermosphere  Ionosphere  Mesosphere  Energetics  and  Dynamics 

TIME-GCM 

Thermosphere  Ionosphere  Mesosphere  Electrodynamic  -  General 
Circulation  Model 

TOC 

Top-of-canopy 

Tk 

Kinetic  temperature 

Tvib 

Vibrational  temperature 
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UARS 

Upper  Atmosphere  Research  Satellite 

UT 

Universal  time 

UVISI 

Ultraviolet  and  Visible  Imagers  and  Spectrographic  Imagers 

VER 

Volume  emission  rate 

VIS 

Visible 

vmr 

Volume  mixing  ratio 

V-V 

Vibration-to- vibration  (transfer) 

WINDII 

Wind  Imaging  Interferometer 
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